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The “Supply” of Physicists 


HROUGH the National Roster of Special- 

ized Personnel we have for the first time 
a detailed breakdown of the professional in- 
terests of practically all of the physicists in the 
United States. Each person, in filling out a 
questionnaire, was asked to indicate in the order 
of his interest or competence five rather narrow 
fields of physics. The sub-divisions from which 
his selections were made were chosen in con- 
sultation with a considerable number of able 
physicists. The charts printed on this and the 
next two pages represent a summary of data 
supplied by the Office of the Roster from all of 
the questionnaires in the field of physics. Since 
it is impossible to evaluate the relative interest 
in the various fields chosen by a given person, 
the charts were drawn merely indicating the 
number of people who have given each field as 
their first interest, second interest and third 
interest. Together, the charts indicate the 
“supply” of physicists skilled in each field. 
Another analysis which 
is now being made by 
the American Institute 
of Physics will indi- 
cate the number of 
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ELECTRICITY & MAGNETISM 


A few words regarding the origin of the 
National Roster of Specialized Personnel and its 
operation are in order. The National Roster is a 
program jointly administered by the National 
Resources Planning Board and the United States 
Civil Service Commission and was formally 
brought into existence one year ago. Dr. Leonard 
Carmichael is director of the program and Mr. 
James C. O’Brien is its executive officer. It was 
formed with two basic principles in mind; first, 
to have a roster listing the qualifications of the 
specialized individual and his availability for 
appointment to the Federal service, and second, 
to conserve the brain power of the nation. — 

The National Roster has been left on a purely 
voluntary basis; there has been no attempt to 
conscript or draft any individual. When there 
are openings in the Federal service where an 
individual registered with the National Roster 
can be of service to the nation, he is offered the 
opportunity of rendering this service, either as a 

consultant or full-time 


eo official. Whether he 
gS seaneiet wishes to accept the 
ME est coe offer or not is entirely 
his own choice; no 
pressure of any type 


is brought to bear on 
any individual. 

As an example of 
this, the following will 
illustrate the methods 
used by the Roster in 
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selecting personnel: 


FIELDS OF SPECIALIZATION 


655 
















































ola & 
ran i 
” H [] 
Vv i 
oO} | | 
i i ° | 
l | : 
uw 1man $ 
re i asl [ i 
a : Nf 3| v 
WJ 2 : Ni . 
ul 2: SBE 3 
S 82 43H 5 2 
5 cme8-4 : | 
y =" - . . 
é 5 fe] 
MECHANICS HEAT 

J HL 
_ | 4 
Ld ie N - 
al | Big | 
>| NNR OF : 
r SN a || 3 
Won SES Ly © ne 
be 3} ‘ s Ele N in HE 
ols) BA oo WA. Ui 
als foot MADE CE Hs 

s -— tt OM fy AN od Fe a3 
=I a he - Spe 
2|: am: aint fT BS Bc bie 

2 Ba Sd eg i ah) < > 

MOLECULAR, ATOMIC 
SOUND & NUCLEAR’ PHYSICS 


(Because of the confidential nature of numerous 
appointments, it is necessary to use a hypotheti- 
cal case.) 

The Navy Department requests the National 
Roster to submit the names of four physicists 
skilled in geometrical and wave optics for in- 
strument design. The request is forwarded by 
the National Roster to its Machine Unit for 
a preliminary selection. The skills of all indi- 
viduals registered with the National Roster 
have been coded and punched on tabu- 


lating cards, in order to speed up the | 


selection process by the use of machines 
to the point where a qualitative analy- 
sis is necessary. The machines then 
select all cards coded 33.3.20, which 
is the code for physicists skilled in 
geometrical and wave optics. The names 
of these persons are sent to the File 





& MAGNETIC MEASUREMENTS 


Unit, where the questionnaires sub- 
mitted by the individuals are removed 
from the file and sent to the Certifica- 
tion Unit for qualitative review. On the 
supposition that the Machine Unit has 


NUMBER OF PHYSICISTS 


ELECTRICAL 
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selected a group of two hundred names, 
the Certification Unit, after an analyti- 
cal review, may select the question- 
naires of thirty-five individuals who meet 
the requirements of the Navy Depart- 
ment’s request. There may be numerous 
reasons why the remaining one hundred 
and sixty-five names were rejected for 
this particular selection. One of the 
major reasons may be that the individ- 
ual is employed either by private in- 
dustry or another Federal agency, where 
he would be of more value to the defense 
effort than in the vacant position in the 
Navy Department. 

The thirty-five remaining names are 
certified to the Navy Department, in 
order that it may contact these persons 
concerning their availability for ap- 
pointment to one of the four positions 
open. The foregoing procedures require, 
generally, less than twenty-four hours, 
from the time the request is received 
by the National Roster until the time 
the names are certified to the request- 
ing agency. 

In the science field the procedure of certifying 
names has been improved by the formation under 
a contract with N.D.R.C. of an Office of Scientific 
Personnel by the National Research Council. 
Dr. Henry A. Barton, Director of the American 
Institute of Physics, is chairman of the govern- 
ing committee of this organization, while Dr. 
Joseph C. Morris is its director. The Office of 
Scientific Personnel is finding that the majority 
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of the calls it receives are for physicists 
and those in related fields in which 
there is such a great shortage of trained 
men. 

As a sidelight, it is interesting to note 
that the first request for personnel re- 
ceived by the National Roster was for 
physicists and the last request received 
as the fiscal year closed June 30, 1941, 
was also from this field. Up until June 
30, 1941, the National Roster had sub- 
mitted the names of more than six 
hundred physicists to various Federal 
and affiliated agencies for consideration 
for appointment. The same demand 
is being made on the Roster for individuals 
in other fields of specialization. It is urged 
that not only physicists but all individuals 
with scientific or professional skills write to the 
National Roster of Scientific and Specialized 
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BORDERLINE FIELDS 


Personnel, Washington, D. C. for questionnaires, 
in order that they may be registered in this 
census of the “brains of the Nation.’’ When 
writing, state what your field of specialization 
is, that is, physicist, chemist, radio engineer, 
etc., so that the appropriate ‘‘check list’’ may be 
sent with the questionnaire. If you have scien- 
tific training and are interested in doing defense 
work, write to the Office of Scientific Personnel, 
National Research Council, 2101 Constitution 
Avenue, Washington, D. C. so that you may be 
placed on the list of available personnel. An 
individual registered in the Roster may never 
be called upon to serve, but his own particular 
skill may be the one necessary to strengthen 
our defense program and thus protect the 
Nation. 
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Résumés of Recent Research 








A Theory of the 
Melting Process 


The most obvious 
properties which distin- 
guish a crystalline solid 
from a liquid are its ability to support a shearing 
stress and the ordered arrangement of its mole- 
cules, demonstrated by x-ray diffraction. In the 
fusion process these properties are destroyed. In 
1939 Born' suggested the vanishing of the 
rigidity modulus as the criterion of fusion. In the 
same year Lennard-Jones and Devonshire’ pro- 
posed a theory of fusion based upon the idea of an 
order-disorder transformation normal 
and interstitial sites of a crystal lattice. Recently, 
John G. Kirkwood and Elizabeth Monroe at 
Cornell University have presented a theory of 
fusion, based upon the order-disorder criterion, 
in which the molecular configuration is not 
artificially restricted to the sites of a basic 
lattice. 


between 


Using classical statistical mechanics, Kirkwood 
and Monroe formulate an integral equation for 
the distribution function of average molecular 
density in a region occupied by a system of 
molecules. They investigate periodic solutions of 
the integral equation subject to the restraint of 
macroscopic homogeneity. When the distribution 
function with the period of a specified space 
lattice is expanded in a Fourier series, the authors 
find that outside a certain temperature-volume 
domain the series reduces to its constant term in 
which all characteristics of the lattice have 
disappeared. The corresponding thermodynamic 
system no longer possesses long range crystalline 
order and is, therefore, identified with the liquid 
phase. 

With the aid of certain approximations the 
authors obtain explicit solutions of their equation 
for the face-centered cubic lattice. They find that 
there should be no critical temperature for the 
liquid-solid transition for this lattice. With the 
Lennard-Jones potential of interatomic force, 
they make numerical calculations of the molar 
fusion parameters of argon. The calculated 
entropy of fusion and volume increment of fusion 
at 83.9°K and one atmosphere are 3.48 E.U. and 
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3.25 cc. Their values compare favorably with the 
observed fusion parameters, 3.36 E.U.and 3.53 cc. 
1M. Born, J. Chem. Phys. 7, 591 (1939). 


? Lennard-Jones and Devonshire, Proc. Roy. Soc. A170, 
464 (1939). 


3 J. G. Kirkwood and E. Monroe, J. Chem. Phys. 9, 514 
(1941). 


Reflection of Light 
from Ceramic Tiles 


In the design of light- 
ing systems for rooms, 
the illuminating engineer 
should have data on the light reflection properties 
of the materials to be used for floor, walls, ceiling, 
and furniture. The reflection factors and colors of 
these materials have a marked effect on the 
amount and spectral quality of the light reaching 
the work. What is even more important, they 
play a dominant role in fixing the brightness 
distribution throughout the visual field, which 
recent research has shown to have an important 
bearing on visual fatigue. 

Hitherto, no comprehensive data have been 
available on the amount and color of light 
reflected from commercially available materials, 
and it is generally impracticable for the engineer 
to make tests on all the materials to be used in a 
given room or building. To meet this need, a 
large number of representative materials have 
been tested at the Massachusetts Institute of 
Technology, and data from these tests have been 
computed and tabulated. The results consist of 
curves of reflection factor vs. wave-length, ob- 
tained automatically on the Hardy recording 
spectrophotometer, and tables of total reflection 
factors and trichromatic coefficients for the 
samples when illuminated with light from incan- 
descent lamps and with daylight (C.I.E. Stand- 
ard Illuminants A and C). 

A recent paper! gives such results for 91 samples 
of ceramic tiles, while forthcoming papers will 
give data on other materials. The tables and 
graphs should be of value to lighting experts and 
interior decorators. Also, since ceramic tiles are 
very permanent in reflecting properties, and are 
carefully controlled for uniformity, the tabulation 
allows tile samples to be used as rough comparison 
standards for judging the color and reflection 
factors of other materials. 


1 Parry Moon, ‘‘Colors of ceramic tiles,’’ J. Opt. Soc. 
Am. 31, 482 (1941). 
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Tri-University Operation of McDonald Observatory 


On July 12 the University of Indiana joined the Univer- 
sities of Chicago and Texas in their operation of the 
McDonald Observatory, located on Mt. Locke in West 
Texas. This observatory is the second largest in the world, 
its 82-inch mirror being surpassed only by the 100-inch 
reflector at Mt. Wilson. 

* 


Summer Program at Brown University 


A program entitled “‘Advanced Instruction and Research 
in Mechanics” has been in progress during the summer at 
Brown University, Providence, Rhode Island. Among those 
giving special lectures were the following: 

PROFESSOR NORBERT WIENER, Massachusetts Institute of Technology. 

Dr. A. NADAI, Westinghouse Electric and Manufacturing Company. 

Dr. HILLEL Poritsky, General Electric Company. 

PROFESSOR S. TIMOSHENKO, Stanford University. 

PRoFEsSOR J. N. Goopier, Cornell University, Buckling and Vibra- 
tion of Bars Involving Simultaneous Bending and Torsion. 

Mr. R. M. Foster, Bell Telephone Laboratories, Electrical Networks 
and their Analog in the General Dynamical Problem. 

PROFESSOR R. D. Courant, New York University, The Rayleigh- 
Ritz Method. 

DR. THEODORE THEODORSEN, Chief, Physical Research Division, 
National Advisory Committee for Aeronautics, Langley Field, Intro- 
duction to the Theory of Wing Flutters. 

Dr. THORNTON C. Fry, Mathematical Research Director, Bell Tele- 
phone Laboratories, Mathematical Research in the Industries. 

PROFESSOR R. B. LINDSAY, Brown University, Filtration of Elastic 
Waves in Fluids and Solids. 


* 
Additions to Staff at Armour Research Foundation 


The appointment of Dr. Raymond G. Spencer and Dr. 
Carl C. Gamertsfelder to the staff of the Light and Optics 
division of the Armour Research Foundation has been 
announced. 

Dr. Spencer, who is conducting x-ray diffraction research, 
is on leave of absence from Albion College, Albion, 
Michigan, where he is Professor of Physics and head of the 
Physics Department. Dr. Gamertsfelder has completed his 
doctorate at the University of Missouri and is working in 
spectroscopy and x-ray investigations. 


* 
Meeting of Scientific Photographers 


The Biological Photographic Association, an_ inter- 
national group of photographers in the natural sciences, 
will hold its eleventh annual meeting in the Hotel Buffalo, 
Buffalo, New York, September 11-13, 1941. 

The program emphasizes the latest methods and proc- 
esses which are likely to contribute to the national defense. 
The annual salon will again be a prominent feature of the 
meeting. Visitors are welcome. For further information 
write the Secretary of the Biological Photographic Associ- 
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ation, University Office, Magee Hospital, Pittsburgh, 
Pennsylvania. 


* 4 
Defense Courses in Physics 


The nation-wide plan for emergency defense courses in 
engineering at the college level, supported by government 
funds and administered by the U. S. Office of Education, 
has been expanded in scope for the fiscal vear which began 
July 1. It now provides for courses in physics, chemistry, 
and production management, as well as the several 
branches of engineering included in the program last year. 
The Advisory Committee of the Office which is composed of 
prominent educators has likewise been enlarged. Dean 
Homer L. Dodge of the University of Oklahoma has been 
appointed to the Committee to represent the field of 
physics in its deliberations. Another member of the 
Committee whose interest lies largely in the field of 
physics is Professor F. L. Bishop of the University of 
Pittsburgh. The Chairman of the Committee is Dean A. A. 
Potter of Purdue University. Inquiries concerning the 
program should be addressed to Commissioner John W. 
Studebaker of the Office of Education, Federal Security 
Agency, Washington, D. C. 


* 


Conference on Visual Problems 


A conference on visual problems was held at Ohio State 
University, June 18-20. The following took part in the 
conference: 


Sevic Hecut, Columbia University, Energy and Vision. 

H. K. HartTLine, University of Pennsylvania, The Discharge of 
Nerve Impulses by the Visual Sense Cells under Conditions of Light 
and Dark Adaptation. 

CLARENCE H. GranaM, Brown University, The Influence of Adapting 
Intensity on the Excitability Cycle of the Photoreceptor Unit. 

S. Howarp BarTLeEy, Washington University, The Optic Nerve 
Discharge and Its Relation to Other Neurophysiological Records in the 
Study of Vision. 

WALTER R. MILEs, Yale University, The Steady Potential of the 
Human Eye. 

Eart M. Lowry AND Davip MacApam, Eastman Kodak Company, 
Visual Sensitometry. 

WOLFGANG KOHLER, Swarthmore College, Fatigue After-Effects in 
Vision. 

J. P. NAFE, Washington University, Adaptation. 

FRANK A. GELDARD, University of Virginia, Persistence of Vision. 

StipNEY M. NEWHALL, Johns Hopkins University, Simultaneous 
Contrast. 

Harry HEtson, Bryn Mawr College, Interrelations of Hue, Lightness 
and Saturation as a Function of Conditions and Modes of Viewing with 
Implications for Measurement and Specifications of Colors. 

DEANE B. Jupp, National Bureau of Standards, Current Research 
Bearing on Choice of a Standard Color System. 

F. L. Drwmick, Hobart College, A Test of Color Aptitude. 

W. B. Lancaster, Dartmouth Eye Institute, Anomalous Spatial 
Localization Associated with Phorias. 

MATTHEW LUCKIESH AND FRANK K. Moss, Nela Research Labo- 
ratories, A New Method of Subjective Refraction Involving Identical 
Techniques in Static and Dynamic Tests. 

FRANK K. Moss AND MATTHEW LuUCKIESH, Nela Research Labo- 
ratories, The Rate of Involuntary Blinking as a Criterion of Relative 
Ease of Seeing. 

FRANCES ILG, Yale University, Directions and Eye Movements in 
Young Children. 

A. M. SKEFFINGTON, Round Table on Problems of Refraction. 
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Ignitor Characteristics 


E. G. F. ARNotTT 
Research Department, Westinghouse Lamp Division, Bloomfield, New Jersey 
(Received December 5, 1940) 


This article reports an investigation of the current and 
voltage characteristics of ‘‘ignitor’’ electrodes in a mercury 
pool rectifier (ignitron). Experiments were made on square 
rods of smooth homogeneous materials of various sizes 
and resistivities with the following results: (1) A frequency 
distribution curve is found for the voltage required by the 
ignitor which is approximately a Normal or Gaussian dis- 
tribution. (2) The total current required for the ignition of 
an arc, divided by the perimeter is a constant, for a given 
material. (3) The total current required per unit length of 


INTRODUCTION 


OR a good many years, the problem of the 
initiation of an arc has been of considerable 
interest to the electrical industry. Various 
methods have been used, such as the thermionic 
cathode, separating contacts, high frequency 
discharge, etc. All of these methods were 
unsatisfactory for the recurrent rapid production 
of high current arcs, where only low voltage 
was available. 

About ten years ago, Slepian and Ludwig! 
discovered that an arc could be initiated by 
passing a current through a rod of high resistance 
material dipping into a mercury pool. 

The current and voltage required for the 
initiation of the arc were not excessive, being of 
the order of 10 amperes and 100 volts. This arc 
starter was developed to a point where it could 
be used commercially and is at present in use in 
various types of tubes manufactured for welding 
and power rectifier applications. Tubes using 
this method for initiating the arc are called 
ignitrons and the starting electrode itself is 
called an ignitor. 

Figure 1 shows a section of a typical sealed-off 
metal ignitron capable of delivering 140 amperes 
average current and several thousand amperes 


!A. I. E. E. Winter Convention, January 23-27, 1932. 
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perimeter times the square root of the resistivity is a 
constant. (4) For identical ignitors operated in parallel, 
the total current required is equal to the square root of the 
number of ignitors times the current required for one 
ignitor. (5) An expression is derived for the fraction of the 
total current which is useful in starting the arc and this is 
found to be inversely proportional to the square root of the 
width. This derivation is verified by the behavior of ignitors 
operated in parallel. 


maximum R.M.S. current. The ignitor is shown 
dipping about }” into the mercury pool cathode. 
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Fic. 1. Section of an ignitron tube. 
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Fic. 2. Test circuit. 


Due to the agitation of the surface of the 
mercury pool, there is considerable variation in 
the current and voltage required by the ignitor. 
It is the purpose of the present paper to discuss 
these current and voltage characteristics and 
to show what relations exist between the 
characteristics and the size and resistivity of 
the ignitor material. 


METHOD OF TEST 


The test circuit is given in Fig. 2 which shows 
schematically the ignitor mounted in a test 
bottle which has a mercury pool and a graphite 
anode. The ignitor is connected to the anode 
supply through a 4-ohm resistance and a hot 
cathode rectifier tube, or a thyratron, so that 
current flows through the ignitor to the mercury 
pool only when the ignitor and anode are both 
positive with respect to the pool. The peak 
voltmeter shown can be connected either across 
the ignitor itself, to measure the ignitor peak 
voltage, or across the 4-ohm resistance to 
measure the ignitor peak current. Figure 3 shows 


Fic. 3. Oscillograms showing 
peak currents (top trace) and 
peak voltages (middle trace) of 
an ignitor and starting voltage 
(bottom trace) of the ignitron. 


ee, 





an oscillogram taken which records the ignitor 
voltage and current, and the tube voltage drop 
for about seventy consecutive cycles. This figure 
indicates the variable nature of the quantities 
which are to be measured. 

Since the temperature of the mercury pool and 
also of the tube walls plays an important part 
in determining the voltage and current required 
by the ignitor, the test tube was operated at an 
average anode current of about 20 amperes for 
at least five minutes before each test reading 
was made. Two methods of measuring the 
characteristics have been used; the first is to 
take oscillograms of the ignitor voltage and 
current, and then to measure the values from 
the films. This method has been extremely 
useful, though somewhat tedious, in determining 
the voltage and current distribution for a single 
ignitor. 

The second method, however, is one which is 
more practical, though less definite, where a 
large number of ignitors are to be tested. This 
consists of using the peak voltmeter and adjust- 
ing it so that it records a voltage or current 
above the setting, less than once in 30 seconds. 
For example, if the peak voltmeter is set at 100 
volts and the ignitor does not require more than 
this value in a period of 30 seconds, then the 
probability of this ignitor requiring more than 
the 100 volts to fire is less than 1 in 1800. This 
method is used in testing ignitors to determine 
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whether the characteristics are suitable for use 


in tubes. It has been used in the following work 
where the determined for 
various ignitors of different sizes and materials. 


characteristics are 


INDIVIDUAL CHARACTERISTICS 


Returning to the ignitor current and voltage 
traces shown in Fig. 2, it can be seen that, in 
general, on the cycles where a high current is 
required to fire, a high voltage is also necessary, 
and for a low current only a low voltage is 
required. This suggests that the resistance of 
the ignitor remains essentially constant from 
one cycle to the next. This is indeed found to 
be the case where the agitation of the pool is 
not too violent. 

A large number of oscillograms were taken 
under various conditions. About 8000 cycles 
were measured visually. The results from several 
films, which were taken as nearly as possible 
under equilibrium conditions, after operating for 
more than an hour, were combined, and a 
voltage distribution curve was plotted. 656 
cycles were used and the data were handled 
statistically to determine the type of theoretical 
curve which would give the best fit to the 
experimental points. 

Figure 4 shows the theoretical curves and the 
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Fic. 4. (1) Voltage distribution curve for an ignitor. 
(2) Area under (1) from — « to V. (3) Calculated distribu- 
tion curve for 2 ignitors in parallel. (4) Calculated distri- 
bution curve for 4 ignitors in parallel. 


experimental points. Curve 1 is the voltage 
distribution curve and shows the fraction of the 
total number of times the ignitor fired in each 
voltage interval of 5 volts. It is a plot of the 
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equation 
1 
a oe 
Xexp [—(V—54.91)?/28.72? AV, (1) 


where P(V) is the probability of firing in the 
interval AV. 


TABLE I. 

V | Ors. f Vy V2 Vsf Vif 

25 0 0 0 0 0 
30 1 30 900 27,000 810,000 
35 12 420 14,700 514,500 18,007,500 
40 38 1,520 60,800 2,432,000 97,280,000 
45 76 | 3,420) 153,900 6,925,500 311,647,500 
50 133 6,650 | 332,500 16,625,000 831,250,009 
55 152 8,360 459,800 25,289,000 | 1,390,895,000 
60 118 7,080 | 424,800 25,488,000 | 1,529,280,000 
65 69 | 4,485); 291,525 18,949,125 | 1,231,693,125 
70 44 3,080 215,600 15,092,000 | 1,056,440,000 
75 13 975 73,125 5,484,376 411,328,125 
80 0 0 0 0 0 


656 36,020 |2,027,650 


116,826,500 | 6,878,631,250 


(8,)!=0.040+0.096, 


1 ’ = Jmin=9.9, 
2+0.24, Bs=2.68 +0.19, 


Juax 2.1. 


Curve 2 is the area under Curve 1 and shows 
the fraction of the total number of times the 
ignitor fired below a voltage V. It is given by 


1 4 —(V—354.91)* 
P< Vj}e—— exp| marenewaeen - av. 
a(2m)' J_, 28.72? (2) 


The other curves, 3 and 4, will be discussed 
later. 

Table I shows the frequency of occurrence in 
the different voltage intervals and how these 
frequencies are used in determining the average 
(V), the standard deviation (¢) and the other 
quantities necessary for determining the type 
of distribution represented by the experimental 
results.2 For a true Gaussian distribution (8;)! 
and J should both be zero. The values obtained 
show that a fit might be had with either a 
Poisson distribution or the binominal law. Since 
the experimental curve is very nearly sym- 
metrical and the Gaussian distribution is at 
times a good approximation to either of these 
two distributions, it was tried first with the 


27. C. Fry, Probability and its Engineering Uses (1928), 
Chapter 9. 
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results shown in Table II. This gives the devia- 
tions of the experimental frequencies from those 
calculated from a Gaussian distribution deter- 
mined from the value of V and o. From these 
deviations, the goodness of fit can be determined, 
and is given by P(> x’) which gives the proba- 
bility of a repetition of the experiment showing 
values further from the theoretical distribution 
than the experimental values already obtained. 
The value of P( > x?) was found to be .55 showing 
a much better fit than would normally be 
expected from such data. In plotting the data, 
it was found that the observer had a personal 
preference for numbers ending in 2, 3, 7 and 8. 
This was eliminated by selecting the voltage 
intervals between numbers ending in 2.5 and 
7.5 volts. 

A similar type of distribution curve was found 
for the ignitor current though sufficient data 
were not obtained to give a complete analysis. 

The reason for the existence of the voltage 
and current distribution curves is not thoroughly 
understood. One suggested explanation is that 
it is due to a general change in mercury level due 


to the agitation of the pool. That this explanation 
is not correct can be seen from the oscillograms 
of current and voltage where it is shown that 
the high voltages occur with the high currents. 
If an ignitor, which has been tested at a given 
immersion, is immersed further into the pool, 
it is found that while the current is increased 
the voltage required is decreased. This clearly 
shows that the deviation from the average 
depends only slightly, if at all, on the general 
mercury level for low anode currents. At high 
anode currents, where the agitation is increased, 
there is some evidence of both factors occurring 
simultaneously. 

The author’s theory is that the agitation of the 
pool causes random fluctuations in the shape of 
the mercury meniscus, which affect the formation 
of the required voltage gradient or energy 
density. This theory is borne out by the results 
obtained by W. E. Berkey on an ignitor immersed 
in a solid tin cathode. In this case, both the 
current and voltage deviations from the average 
were extremely small. Fortunately, it seems 
probable that a recording counting circuit can 




















TABLE II. 
Yo w 
pPu(V) -——|J exp (y?2/2)dy — | ; exp (—y? nay]. exp ( —y?/2) =exp [—(V1 —54.91)2/2¢?] 
y ae | 
; : mae Vi—54.91 ' : ” m | 
J (V1) (Vi —54.91) ( a ) mon” i vied (—y?/2)dy p(V) Exp.f | Oss. f | DEv. | Div 
30 0.005 a 1 —2.3 1.600 
32.55 | —22.41 —2.57 0.005 
35 .018 11.8 12 | +0.2 .004 
37.5 —17.41 — 2.00 .023 | 
40 055 36.1 38 | +1.9 .100 
42.5 — 12.41 — 1.42 .078 
45 | 122 80.0 | 76 —4.0 .200 
| 47.5 — 741 —0.83 .200 
50 197 | 129.2 133 +3.8 | 115 
| §25 — 2.41 —0.28 397 | 
55 223 | 146.0 | 152 | +6.0 | .247 
57.5 + 2.59 +0.30 .620 
60 .190 124.5 | 118 | -—65 | .338 
62.5 + 7.59 +0.87 .810 
65 115 | 75.5 69 —6.5 | .560 
67.5 +12.59 +1.44 .925 | 
70 .053 34.8 44 +9.2 | 2.430 
72.5 +17.59 +2.02 .978 | 
7S | 017 | 11.2)} = 13 
|. Seon | +22.59 +2.59 .995 —1.5 155 
80. | | .005 3.3} | 0 
| 20 00 | a 1.000 
| 1.000 | | 5.750" 
| 
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| S'=7, P(> x*)=0.55 
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be made which in the future, will give in a few 
minutes, the data which required many hours of 
measuring to obtain before. 


GENERAL CHARACTERISTICS 


In the earlier work on ignitors, current and 
voltage readings were recorded for different 
materials. It was noticed that, while there was 
a great deal of variation in the results, most of 
the values obtained showed a definite tendency 
for the current to increase as the voltage de- 
creased. When the voltage versus current was 
plotted, it was found to approach a rectangular 
hyperbola of the form VJ=K. Figure 5 shows 
the approximate range of the experimental 
points and the curve which seems to represent 
the average. The values of current and voltage 
found in these measurements were all taken 
with a peak voltmeter in the manner already 
described. 

Straight rods of square cross section accurately 
ground to within +.001” were then obtained of 
different resistivities and having smooth surfaces. 
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Fic. 5. Approximate range of voltage and current char- 


acteristics of a large number of ignitors showing the wide 
range of values. 
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TABLE IIT. | 
MaTL.| p | 
No. |@/In3| d V I lo pe! | Top! |\V/p'iV(4a)3 
3 | 1.73 |1/16’") 198} 8.3] 33.2 |1.32 143.7 |150 | 
11 | 2.27 |1/16’"| 204| 7.2) 28.8 |1.51 |43.5 |135 | 
5 | 2.25 |1/16’"| 206 | 8.0 | 32.0 |1.5 |48.0 |137 | 
5 | 2.45 |1/16’"| 207); 7.4) 29.5 1.56 46.0 132 | 
| | | 
1 1/16’’|*129 |*15.4| 61.5 |0.82 |50.5 157 | 64.5 
1 | 0.64 |5/64’") 114} 19.5) 62.5 |0.80 |50.0 |143 | 65.0 
1 | 0.71 |3/32’"| 107 | 22.7) 60.5 |0.84 |50.8 |124 | 65.5 
1 | 0.69 |7/64’| 97! 26.8) 61.5 |0.83 |51.0|117| 64.0 
| | 
2 | 0.21 |1/32'"| 74) 14.0/|112.0 |0.46 51.5 | 
2 | 0.22 |1/16’") 70) 27.8 111.0 |0.47 |52.2 |149 | 


* Calculated. 


All the samples used on these tests were 1}” 
long and were immersed }” in the mercury 
with about ?” between the mercury and the 
ignitor holder. On testing samples of different 
sizes, it was found that a definite relation existed 
between the currents and the width of the side 
of the rod. This relation was that the current 
divided by the perimeter (or by the length of 
the side) was a constant for a given material as 
is shown in Table ITI. 


I/4d =constant = I». (3) 


This suggests that for a given material, a 
definite linear current density or current per 
inch of perimeter is required to start the arc. 
Each of the values given are the average of those 
obtained for several ignitors of the same size 
and material. 

Three different materials were used in these 
tests and when the voltage was plotted against 
the current per inch, the wide variation similar 
to that which existed in the older current- 
voltage curves, largely disappeared. Figure 6 
shows the curve obtained. It will be noticed 
that the points for the material with p=.68 
ohm/cm* are spread vertically. For this set, it 
was found that the voltage times the square 
root of the width was approximately a constant. 
However, there is not sufficient variation in size 
to make this very conclusive except over a small 
range. The point at 129 volts and 61.5 amperes 
per inch was calculated for this constant for an 
ignitor with a ;,”’ side. This allowed the heavy 
curve to be drawn through three points all for 


is’ ignitors. 
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Fic. 6. Voltage plotted against current per inch for ignitors 
of different sizes and materials. 


The light lines running from the origin to the 
maximum peak voltage and current points are 
drawn to show where the lower firing values due 
to the distribution curve will fall for each 
ignitor. This curve has a different shape from 
that of the earlier one, but when account is 
taken of the fact that the lower resistivity 
material operates at a considerably higher 
temperature, it is believed it will return to a 
curve represented by 


(V—a)(IJy—6) =constant, (4) 


where b is probably =0. 

When it was found that the current per inch 
was constant for different materials, an attempt 
was made to find out what relation existed 
between this current and the resistivity. This is 
complicated by the fact that all the materials 
used have a large negative temperature coeffi- 
cient of resistivity. When an ignitor is operating, 
its temperature may reach several hundred 
degrees centigrade, and is not uniform along its 
length. Thus it is very difficult to determine 
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the hot resistivity near the surface of the 
mercury with any degree of accuracy. However, 
the fact that for different size ignitors of the 
same material J/4d is a constant indicates that 
the temperature of the ignitor very close to the 
mercury surface is determined largely by the 
temperature of the pool itself. When the values 
of cold resistivity were compared with the 
current, it was found that approximately the 
current was inversely proportional to the square 
root of the cold resistivity. This suggests again, 
as the early current voltage curves had, that a 
definite amount of energy is required to form 
the arc. 

Figure 7 shows 1/p’ plotted against J». The 
points are reasonably close to a straight line. 
Since this line should go through the origin, 
the two upper sets of points are progressively 
further from it. It is felt that this is due to 
the resistance versus temperature curves being 
slightly different for the three different materials 
used. 

It should be emphasized that the current and 
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Fic. 7. 1/p' versus the maximum peak current per inch. 
The dotted curve is the apparent straight line which fits 
the points. The solid curve is the straight line which is 
probably more correct if account is taken of the difference 
in the temperature coefficients of the materials. 
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voltage readings being discussed are the maxi- 
mum peak values which are not exceeded in a 
time of 30 seconds. While the above results are 
of considerable practical value in designing 
ignitors, they must be regarded with caution in 
trying to determine the fundamental processes 
which are responsible for the starting of an arc 
by this method. For this purpose, it is felt that 
the average or most probable values are the ones 
which should be used in investigating the 
fundamental processes. This fact must be re- 
membered, since it is not necessarily true that 
the same relations will hold for the most probable 
values given by the distribution curves as hold 
for the maximum peak values. 


IGNITORS OPERATED IN PARALLEL 


Since the firing of an ignitor is of a statistical 
nature and the measurement of voltage and 
current can be interpreted as giving the proba- 
bility of firing, this combined with the relation 
I /4d=K suggests that further information might 
be gained by operating two or more ignitors in 
parallel. This was done with the square ignitors, 
and it was found that the maximum peak current 
required for nearly identical ignitors was given 
approximately by 


T,=n'I,, (5) 


where n is the number of ignitors and J, is the 
current required for one ignitor operated alone. 
From Eq. (3) 
I /4d=K, 


we would expect to get 
[,=nl; (6) 


for nm ignitors in parallel and the voltage should 
remain about the same. 

The following argument, however, shows that 
the observed results are to be expected, at least 
qualitatively. 

Let P=the probability of one ignitor firing 
below V volts. Then (1—P)=the probability 
of one ignitor firing above V volts. If two 
identical ignitors are operated in parallel, then 
(1—P)(1—P) =(1—P)?=the probability of both 
ignitors firing above V volts and P2=1—(1—FP)? 
=the probability of one or the other, or both 
ignitors firing below V volts, P will, in general, 
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be given by 


P,=1-(1-—P)", (7) 


where m is the number of identical ignitors 
operated in parallel. Returning to the case of 
just two ignitors in parallel 


P,.=1-—(1-—P)?, (8) 
dP,=2(1—P)dP. (9) 


Now a plot of dP versus V is just the distribu- 
tion curve shown previously and dP», versus V 
will be the distribution curve for two ignitors. 
This curve and one for n=4 are shown in Fig. 4, 
curves (3) and (4) and both of them are still 
reasonably symmetrical about the most probable 
value and resemble a Gaussian distribution. 

It is seen from curves (1) and (3) that with 
two ignitors, the probability of firing is increased 
for values of V below the most probable value 
for one ignitor, and is decreased above this value 
of V. A similar relation is found between curves 
(3) and (4). 

This analysis is, of course, only justified if 
each ignition is completely independent of the 
previous one. This in general seems to be the 
case since, on a rough test, the values of V are 
found to be distributed at random with respect 
to time. 


THEORY 


The theories proposed for the starting of an 
arc by an ignitor depend on either the formation 
of a large voltage gradient or a high energy 
density at or very near the point where the 
surface of the mercury joins the ignitor. Slepian! 
has shown that the current density will be ex- 
tremely high at the point of contact. However, 
the experimental results obtained suggest that 
not all of the current which passes from the 
ignitor into the mercury is useful in starting the 
arc. 

It is possible to derive an expression for the 
ratio of useful to total ignitor current and this 
will be used with the experimental results to 
determine the effect of surface contact area on 
the total current required to start the arc. The 
experimental values of the total current referred 
to are the maximum peak values obtained as 
described above, with a peak voltmeter. As 
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shown by Fig. 2, a current distribution exists 
similar to that shown for voltage and it is for 
this reason, that the value of the firing current 
is affected less than would be expected by 
change in size of the ignitor and by the number 
of ignitors used together. 

Expressions are given by Maxwell’ for the 
equipotentials and lines of current flow in a thin 
sheet of a conductor of such shape that it 
represents a vertical cross section of an ignitor 
dipping into a mercury pool. It is assumed here 
that these expressions can be used as a suffi- 
ciently good approximation for the three-di- 
mensional current flow to give the ratio of the 
current entering the mercury above any point 
to the total current flowing into the mercury. 

In Fig. 8, x is the distance from the point of 
contact with the mercury 2b=d is the length of 
one side of a square ignitor, and / is the length 
of ignitor immersed in the mercury. 

Then if ¢ is the potential function and y the 
current function, Maxwell gives 


b 1 
«=— log —(exp (2¢) +exp (—2¢)+2 cos 2p), (10) 
2 4 


, ' (exp (¢) —exp (— ¢) tan y) 
y=b tan—! ——_ —, 


exp (¢)—exp (— ¢) 





(11) 


The equipotentials are given by setting 
¢=constant and the lines of current flow by 
setting Y=constant. 

To find the amount of current which enters 
the mercury close to the surface we set g¢=0 
for the surface of the ignitor below the mercury. 
This gives 


b 1 
a hae 2y) (12) 


from which 
¥Y=cos™! exp (x/b). 


Since we are concerned only with the distance 
below the mercury contact, we can write 
Y=cos! exp (—x/b), (13) 


where x is now taken as the distance on the 
surface of the ignitor below the mercury. 


3 Maxwell, Electricity and Magnetism, third edition, 
article 93, p. 297. 
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From this equation, the amount of current 
passing into the mercury within a distance x of 
the surface is given by 


¥z2—wWo=cos exp (—x/b) —cos™ e® 
=cos! exp (—x/bd) 


and the proportion of the total current which 
passes into the mercury within a distance x of 
the surface is 

Y2—Wo cos! exp (—x/b) 


Vi-—Wo cos" exp (—1/6) 


cos~! exp (—«/b) 
= for 1>>7b. 
w/2 





Thus if we call z the useful current and J the 
total ignitor current 


2 
=— cos! exp (—x/b). 
T 


(14) 


mm | =. 


If x is taken as .001”, rh=d as .062”, 1 as .25”, 
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Fic. 8. Equipotentials and lines of current flow in an ignitor 
immersed in mercury. 


then 7=0.27. Since the value of x must be 
extremely small, probably much less than .001” 
exp (—x/b)=1—<x/b approx. 


and 
11 x 


similarly for small values of 7/J 


ri 31? 
cos — = 1 ——— approx. 
2 
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and 


es =z 
sl b 
or 
t=I(8x/m*b)! (15) 


or since d, the width of the ignitor is equal to 7b 
i=I1(8x/md)'. (16) 


This equation shows that the proportion of 
useful current decreases as the width of the 
ignitor increases. Let xo be the distance below 
which current entering the mercury cannot aid 
in starting the arc. Then 

t=1(8x9/ad)'=al /d' (17) 
where 

a= (8xo/m)'=constant. 


Since it was found that the current required to 
fire the ignitor 


IT=4dIo, (3) 
then the useful current will be given by 
t=4dI,ya/d' (18) 
or 
t=4alI od). (19) 


This shows that the useful current required 
to start the arc is proportional to the square 
root of the perimeter of the ignitor for square 
ignitors of a given material. 

In the case of n identical ignitors each of 
width d operating in parallel, it is seen from (17) 
that for each ignitor, the ratio of useful to total 
current will be 

1, I,=a (d,)} (20) 
and 
in/I,=a/(d;)}, (21) 


since this ratio will remain constant regardless 
of the number of ignitors used because they are 
all of the same size. 

But from (19) 


tn = 4aIy(nd,)}, (22) 


since 4nd, is now the total perimeter in contact 
with the mercury and 


in(d,)! 4aIo(nd,)*(d,)' 


n 
a a 
or 


IT, =4Iodn'. (23) 
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Since, in this case, d is a constant, then 
IT, =n'J\. (24) 


That is, the total current required to fire a 
number of identical ignitors in parallel is 
proportional to the square root of the number 
of ignitors. 

This result has been approximately verified 
experimentally, and thus justifies the approxi- 
mations made in deriving the expression for the 
useful current. Dividing the useful current by 4d, 
we get the value of useful current per inch 
around the ignitor which it requires to start 
an arc. 


Io =1 4d =4aI,d' 4d=al, d}, (25) 


since J) is a constant, ip cannot be, and the 
reason for this is to be found in the statistical 
nature of the firing process. 

That this result is to be expected qualitatively 
can be seen by considering the surface of the 
ignitor, at the junction with the mercury, as 
composed of a large number of small points, 
any one of which is ready to start an arc under 
the proper conditions. A certain number of 
these will be found ready to start an arc for some 
value of useful ignitor current per inch, to’. As 
the perimeter of the ignitor is increased, the 
number of points available for firing is increased 
so that the chance of there being at least one 
point ready to fire at a given value of current 
will also be increased. This means that the 
chance of firing at low values of i is increased, 
and therefore, a new value 79’ can be found 
which is less than i’, for which the chance of 
firing is the same as it was for 7%’ with the 
smaller ignitor. 

Thus the total useful current 4dio will not 
increase linearly with d but will increase more 
slowly than this. 


SUMMARY 


It would be unwise at the present stage of 
this work to try to draw too definite conclusions 
from the results. It is felt, however, that the 
following general statements are justified for 
square rod ignitors of smooth homogeneous 
material. 

1, A frequency distribution curve can be 
obtained for the current and voltage required to 
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fire an ignitor. This curve is approximately a 
Gaussian error distribution. 

2. The total current required to fire, divided 
by the perimeter, is a constant for a given 


5. For identical ignitors in parallel, the current 
required to fire is equal to the square root of 
the number of ignitors times the current required 
for one ignitor. 


resistivity. I,=n'I,. 
I/4d =I». 6. The effect of the statistical distribution on 
3. The total curveat per unit leneth of the firing current is proportional to the square 


root of the perimeter. 

Work is being continued in order to determine 
from the distribution curves how the most 
probable values depend on the size and resistivity 
of the ignitors. 


perimeter times the square root of the resistivity 
is a constant 


Ipp'=K. 


4. The proportion of useful current required 
to fire the ignitor is inversely proportional to the 
square root of the width 


The writer wishes to express his appreciation 
to Mr. D. E. Marshall for his continued help and 
criticism and to Dr. J. Slepian for his comments 
during the various stages of the work. 


i/I=a/d'. 





The Structure and Dynamics of Liquids* 


W. E. RosEvEARE, R. E. PowELL AND HENRY EyYRING 
Frick Chemical Laboratory, Princeton University, Princeton, New Jersey 


(Received December 7, 1940) 


By coordinating our information on viscosity, diffusion, melting, and other rate and thermo- 
dynamic properties, we arrive at a detailed picture of liquid structure. Thus we find that a 
liquid is best thought of as a solid to which a large number of empty equilibrium positions 
are added. In fact the expansion on melting, as well as the expansion with a rise in temperature, 
arises almost entirely from this introduction of new equilibrium positions, and only to a minor 
extent from lattice expansion. We shall obtain information as to the number, size and energy 





of formation of these empty lattice points. 


I. THEORY OF MELTING 


EVERAL mathematical treatments of phe- 

nomena in the liquid state have taken as a 
model liquid a binary mixture of holes and 
molecules. Thus Eyring was led to postulate the 
presence of holes in liquids in a theory of the 
law of rectilinear diameters of Cailletet and 
Mathias as well as to explain viscous flow.! 
Lennard-Jones and Devonshire* assume that in 
the process of melting a number of new lattice 
positions called 6-sites equal in number to the 
original a-sites appear. The redistribution of the 


* Presented at the 1940 annual meeting of the Society of 
Rheology. 

1H. Eyring, J. Chem. Phys. 4, 283 (1936). 

2 J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. London A169, 317 (1939); ibid. A170, 464 (1939). 
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atoms between the old a-sites and the new £-sites 
involves an increase in entropy, energy, and 
volume, and with a suitable and reasonable 
choice of parameters excellent agreement with 
experiment has been obtained by these authors. 
Cernuschi and Eyring have considered liquids as 
being made up of a binary mixture of holes and 
molecules, the holes having the same size as 
the molecules. A mathematically equivalent 
treatment has been given by Kirkwood who 
treated critical phenomena as order-disorder 
transformations.‘ 

We will now make the assumption that a 
liquid is a random distribution of molecules over 

3F. Cernuschi and H. Eyring, J. Chem. Phys. 7, 547 


(1939). 
‘J. G. Kirkwood, J. Phys. Chem. 43, 97 (1939). 
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TABLE I. 


V—V,. | Votume or Ion nh ASu aSu 
SUBSTANCE V. VOLUME OF Atom| WN CaALe Oss. 
Mercury 0.0259 0.114 1.363) 3.2 e.u. 2.37 e.u 
Lithium .006 .097 0.37 | 1.58 1.83 
Sodium 0271 .170 0.955) 2.68 1.70 
Potassium 0231 .206 0.672) 2.24 1.75 
Rubidium .0284 .245 0.695) 2.27 1.67 
Cesium 0256 .273 0.562) 2.08 1.69 
Benzene .1338* 0.801) 2.46 8.5 
eh .0677 0.406, 1.68 8.5 
Carbon di- 
sulfide .0673 0.404, 1.67 
Argon 12 0.72 | 2.32 3.19 
Nitrogen .09 0.54 | 2.00 2.69 


*From Bridgman's data on the change of melting point with 


pressure. 


occupied and unoccupied equilibrium positions, 
and that the volume and energy associated with 
one of these equilibrium positions is characteristic 
of the material. The fusion of a mole of mole- 
cules, N, results in the introduction of m, of 
these holes, to give an entropy increase of 


Nh 
—n,k ln ———. (1) 
N+n, 


AS=—WNk In 
N+n, 


If the volume increase on melting is due to the 
introduction of such new equilibrium positions, 
we have: 


n,/N=n'(v—v,)/Vs, (2) 


where v, is the volume of the solid, v is the volume 
of the liquid at the melting point, and the volume 
of an empty equilibrium position is taken as 
1/n’ of the volume of a molecule. A numerical 
value of n’ may be found as follows: Batschinski 
has shown that for a large number of liquids the 
fluidity is proportional to the increase in volume: 
or in our terminology, the fluidity is to a large 
extent determined by the number of empty 
equilibrium positions alone.®® From a study of 
the pressure effect on fluidity, Frisch, Kincaid 
and Eyring have shown that the volume of one 
of these holes required for flow is about § the 
volume of a molecule.’ 

Taking n’ as 6, the observed expansion on 
melting permits us to compute by Eq. (2) the 
number of holes in the liquid at the melting 
point and thus permits a rough computation by 

*A. J. Batschinski, Zeits. f. physik. Chemie 84, 643 
EC Bingham, Fluidity and Plasticity (McGraw-Hill, 


7D. Frisch, H. Eyring and J. F. Kincaid, J. App. Phys. 
11, 75 (1940). 
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Eq. (1) of that part of the entropy of melting 
due to the increase in lattice disorder. Table I 
presents the calculated and experimental results. 
Fluidity data has been interpreted as showing 
that in molten metals the ion is the flowing 
species, so for the metals, 7’ is taken as 6/ (volume 
of ion/volume of atom). 

The calculated entropies of fusion are seen to 
cluster about 2 entropy units: this is just the 
observed entropy of fusion for a large number of 
monatomic substances.* Except for the simplest 
substances, melting will involve other entropy 
changes besides that arising from the increase in 
disorder.'* "' The number of holes just equivalent 
to 2 e.u. is m,=0.54N. The results of Table I 
may be summarized by saying: The entropy of 
melting due to lattice disorder will be reproduced 
by the distribution of 1 mole of molecules at 
random upon 1.54 moles of lattice positions. 
Further, the volume change of melting will be 
reproduced if the volume of a hole is about 3 the 
volume of a molecule. The thermodynamic equa- 
tion dp/dT=AS,,/AV, insures that a_ theory 
which gives correctly the entropy of melting and 
expansion of melting will also give the proper 
dependence of the melting point on pressure. 


II. THE PARTITION FUNCTION FOR A LIQUID 


We can now formulate, in terms of holes, a 
simple partition function for a liquid. The parti- 
tion function for the liquid is intermediate 
between that for the solid and that for the gas. 
In a solid the molecules vibrate in a field of the 
type shown in Fig. 1a, while in a liquid some of 
the degrees of freedom are of the type shown in 





(a) 
Fic. 1. 


® R. Ewell and H. Eyring, J. Chem. Phys. 5, 726 (1937). 

®Cf. Taylor, editor, Treatise on Physical Chemistry 
(Van Nostrand, 1931) second edition. 

1 L. Pauling, Phys. Rev. 36, 430 (1930). 

1 J. Hirschfelder, D. Stevenson and H. Eyring, J. Chem. 
Phys. 5, 896 (1937). 
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Fig. 1b, where D is the activation energy for 
passage to an adjacent vacant equilibrium posi- 
tion. The partition function for one degree of 
freedom of a molecule in a potential field of the 
type shown in Fig. la is: 


fi=(1—exp (—hv/kT))— exp (—E,/3RT), (3) 


where v is the characteristic vibration frequency 
of the solid and E, is the potential energy of the 
system when all the molecules are at their 
equilibrium positions. The number of degrees of 
freedom of this type will be three times the 
number of molecules multiplied by the mole 
fraction N,(v) of filled equilibrium positions. The 
total partition function for all degrees of free- 
dom of this type will be: 


F,=[(1—exp (—hv/kT))“ 
Xexp (—E,/3RT) Y"™, (4) 


For a molecule in a potential field of the type 
shown in Fig. 1b, the partition function in that 
degree of freedom will be that for a simple har- 
monic oscillator up to the internal dissociation 
energy D. Above this we have half the partition 
function for the harmonic oscillator plus the 
partition function of a particle in a box of length 
(v,/n’)' as indicated by the adjacent dotted lines. 
The partition function for a single degree of 
freedom of this type is: 


D a) 
fe-| > exp(—nhv/kT)+} > exp(—nhv/kT) 


nhv=0 nhe=D 
(2amkT)*/ v.\3 
gue —(=) exp (—D eT] 
h n’ 


Xexp(—E,—3D/3RT) (5) 


wn 


or 


, k —} exp (—D/kT) 


1 ~exp ( — he/ kT) . 


(2emkT)! 70, \! 
+exp (—D/kT)— ——(-) | 


h n’ 


Xexp (—E,—3D/3RT). (6) 
The number of degrees of freedom of this type 
will be three times the number of molecules 


multiplied by the mole fraction N2(v) =1— Ni(v) 
of holes. The partition function for degrees of 
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freedom of this type for N molecules will be: 


‘ ( —} exp (—D/kT) 
LN 1-exp (—hv/kT) 


(2mmkT)* sv. \! 
+exp (Dyer) —(~) ) 
h n' 


3N No(v) 
Xexp [—(E,—3D) RT| ‘ (7) 


In addition the partition function contains a 
term corresponding to the entropy of lattice 
disorder : 


(N+n,)! 


3=> 


N! n nf 





(8) 


The complete partition function for N molecules 
is then: 


F= F,F.Fs3. (9) 


This partition function becomes at the lower 
limit that for the solid, at the upper limit that 
for the gas. Since the energy of making a hole is 
given as 3 the energy of vaporization, the parti- 
tion function will give the critical point approxi- 
mately correctly.” 

From the well-known equation 


A=-—kT Inf (10) 
and equations of the type 
P=-—(0A/dV)r (11) 


the thermodynamics of liquids may be derived 
when an expression is found for N,(v). Applica- 
tions of this type of a partition function to 
thermodynamic properties will be given else- 
where.!!* 


III. FLurpiry AT CONSTANT VOLUME 


Bridgman’s data at high pressures shows that, 
for non-associated liquids, the temperature coeffi- 
cient of fluidity at constant volume is insig- 
nificant compared to the temperature coefficient 
at constant pressure.'* These results indicate that 


"a J. Walter and H. Eyring, J. Chem. Phys. 9, 393 
(1941). 

2 J. G. Kirkwood, J. Chem. Phys. 7, 909 (1939). 

3 P. W. Bridgman, The Physics of High Pressures 
(Macmillan, 1931). 
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the fluidity of normal liquids under ordinary 
conditions is to a large extent a function of the 
volume alone. 


As we shall see later, diffusion can be best 
interpreted as the passage of a single molecule 
over a potential barrier into a hole. Thus the 
most likely mechanism of flow is the moving of 
single molecules into adjacent holes. 

For a molecule located at a particular point, 
its average velocity of moving forward will be 
equal to the product of the number of equi- 
librium positions, or holes, to which it may flow 
times the average velocity of the molecule 
parallel to the direction of flow times an expo- 
nential term involving the free energy of activa- 
tion for flowing from one equilibrium position to 
another. This latter activation free energy is 
generally small compared with the free energy 
of forming the hole F,. The number of new 
equilibrium positions available may be set equal 
to the increase in volume over that of the solid, 
i.e., v—v, divided by v, the volume increase 
required per new equilibrium position. The aver- 
age velocity of the molecules parallel to the 
direction of flow is (k7/2xm)!. The average 
distance between equilibrium positions in the 
direction of motion is taken as \ while the dis- 





tance between neighboring molecules in the 
TABLE II. 

SUBSTANCE AH, SUBSTANCE AH, 
Mercury 700 Acetone —110 
Ether 700 n-octane 1180 
Benzene 700 methyl alcohol 1340 
Chlorobenzene 640 ethyl alcohol 3030 
Brombenzene 740 n-propylalcohol 3000 
lodobenzene 780 i-propyl alcohol 3640 
Carbon tetrachloride 530 n-butyl alcohol 3000 
Carbon disulfide 530 t-butyl alcohol 4120 
n-pentane 750 t-amyl alcohol 4390 
i-pentane 460 Glycerine 1820 
Hexane 640 
same direction is A» The distance between 


molecules normal to the direction of flow and 
perpendicular to the velocity gradient is )s. 
The shearing force per square centimeter is f. 
The force acting on a single molecule is fds, 
and it acts to lessen the work of passing over the 
potential barrier through a distance \/2, so that 
in the forward direction the height of the barrier 
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is in effect lowered by the amount of work 
frX2xA3A/ 2, while in the backward direction it is 
raised by the same amount. The number of 


.times that a molecule moves in the forward 


direction in a second may then be written as the 
corresponding specific reaction rate: 


v—v,1sf kT \! 
v, A\2aem 


Xexp [(—AH,+ frad3d/2)/kT]. (12) 
For the backward direction : 
v—v, 1 kT \3 
WH) 
v, A\2rm 
For the fluidity we have: 
Av A\(k;—kR,) 
¢=— = ———__, (14) 
fd fm 
h v—v,1/ kT \3 
=— —— -( -) exp (—A//,/kT) 
fX¥1 vn A\2Zaem 
X2 sinh (frAad3A/2kT). (15) 


For ordinary viscous flow fdaA3\«Kkt, so that 
expanding the exponential, retaining only the 
first power terms and letting \2\3;= V! we have: 


9) 


Ch 


r U—U,z 
o=—ai( —) (2rmkT)-bexp (—AH,/kT). (16) 
i 


If we take \/A; equal to (v,/v)! and v,=2v,/n’, 
the fluidity is given as follows: 


g9— 9, 
6=(om'(* )w’(emk)- 


9) 


Us 





Xexp (—AH,/kT). (17) 


Taking n’ as a constant, we have at constant 
volume for two different temperatures: 


1 
AH, =k(A In $+ 3A In r) /a(—). (18) 


Table II gives values of A//, in calories per mole 
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TABLE III. 





MELT- . —AH+ ) o 
inc |® exp ( kT ; CALcu-|  Os- 

SUBSTANCE POINT CALCULATED LATED | SERVED 
Mercury 234°K | 51.5 j11 48.8 
Cadmium 250 | 45 | 66 
Lead | 601 41 | 68 
Tin | 505 55 | 24 
Bismuth 544 44 | 54 
Nitrogen 63 | 484 | 339 
Argon 84 633 } | 357 
Benzene 278 | 252 71 | 121 
Carbon disulfide 160 | 242 46 | 59 
Carbon tetrachloride 250 206 71 49.2 
Chlorobenzene 228 | 378 92 | 44.4 
Bromobenzene 243 227 49 37 
lodobenzene 242 175 30 22 
n-pentane | 142 | 410 \31 24 
n-hexane 178 | _ 356 |58 50.8 
n-heptane 283 302 37.2 
n-octane 216 | 314 24 | 56.2 
n-decane 243 350 | 47.6 
o-xylene 244 328 | 56.4 
Methyl acetate 175 354 31.4 
Chlorine 170 178 110 
Bromine 266 131 69.0 
lodine 386 101 44.8 
Sodium chloride 1077 61.9 
Potassium nitrate 35.7 
Sodium nitrate 32.9 
Ethyl alcohol 156 322 050 | .031 
n-propyl alcohol 146 314 068 | .034 
l-amy]l alcohol 156 320 00003 .00007 
Glycerine 290 314 | .087 .082 


determined from this formula using fluidities at 
30° and 75° at constant volume. These fluidities 
were obtained by graphical interpolation of 
Bridgman’s data in the lower pressure range 
where the temperature coefficient of fluidity at 
constant volume increases with increase in den- 
sity. The uncertainty in the determined values 
of AH, is about 100 calories. The normal sub- 
stances in the first column of Table II have 
values of A/T, constant within about the experi- 
mental error. Acetone like water has a negative 
value of AH7, which may be interpreted as due 
to a change in structure or coordination number 
with temperature at constant volume. Most of 
the alcohols have activation energies at constant 
volume corresponding in magnitude to the energy 
required to break a hydrogen bond. 

At the melting point we may substitute Eq. 
(2) into Eq. (17), set m,=0.54N, and obtain: 


(—AH,) 


dm = 0.54(v,0)'(2amkT)—} exp —" (19) 


Table III gives the observed values of the 
fluidities at the freezing point for various types 
of substances, and those calculated from Eq. (19). 
The calculated values agree with the experi- 
mental values within a factor of about two. The 
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uncertainty in the values of A//, is sufficient to 
account for all divergence between the calculated 
and experimental values except in the case of 
mercury, which is four times more fluid than 
calculated. This discrepancy is removed if in 
Eq. (19) we take n’ as 6/(volume of ion/volume 
of atom) for metals as we should if the ions 
rather than the atoms are the units of flow. 
The value of the exponential factor for most 
substances varies from four to six, while for the 
hydrogen bonded substances the value is very 
much larger. It is 10’ for isoamyl alcohol. The 
table also indicates that except for the sub- 
stances melting below 100°K the fluidities at the 
freezing point tend toward a value of about fifty 
for normal substances. 

If v,, n’ and AH, are independent of volume, 
then @7'/[v! exp (—AH,/kT) | should be a linear 
function of the volume, according to Eq. (17). 
Figure 2 shows this to be true for benzene, 
carbon tetrachloride, carbon disulfide, chloro- 
benzene, isopentane and n-butyl alcohol. The 
range of data plotted for benzene is 10—185°C, 


for carbon tetrachloride 0—-65°C, for carbon 
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disulfide 0-46°C, for chlorobenzene 0—120°C, 
for isopentane 0-30°C and for n-butyl alcohol 
0-38°C. In Fig. 3, mercury gives a straight line 
from —20 to 100°C while showing very marked 
curvature from 100 to 350°. In this latter range 
it is probable that the value of n’ is decreasing 
toward the value for normal nonmetallic liquids. 
In Fig. 3, ether, ethyl alcohol and glycerine do 
not give straight lines, suggesting a change in 
the liquid structure. (In Figs. 2 and 3 the scales 
of coordinates are for benzene and mercury 
respectively, while the values for the other 
substances are in arbitrary units and with the 
axes displaced to fit the graph.) 


IV. FLurpiry AT CONSTANT PRESSURE 


Since most of the available fluidity data are at 
constant pressure, instead of at constant volume, 
it will be useful to see what results the hole 
theory of liquid structure yields in this case. 

The equation for fluidity at constant pressure 
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may be written® 


1 v 
-=o=— exp (—AF*/RT), (20) 
n Nh 

=— exp (AS*/R) (exp —AH*/RT), (21) 


Nh 


where the term AH now contains both the heat 
of activation for forming the hole and the heat 
of activation at constant volume, and the term 
AS* includes any other factors. This is, of course, 
only an alternative way of writing previous 
expressions. 

The free energy of activation AF* may be 
estimated as some fraction of the total free 
energy change AF of an analogous reaction. In 
order to make a hole, or new equilibrium position, 
the same kind of bonds must be broken as are 
broken in vaporization; however, the work RT 
per mole of expansion against the atmosphere is 
not needed, neither is the entropy of expansion 
into the vapor state available. Thus the thermo- 
dynamic quantity with which AF* may be corre- 
lated is 


AF yap t+ TASvap — RT = AA vay — RT = AE yay. 


Ewell and Eyring have pointed out that AH* 
is about 4AE,,, for a number of liquids. Fig. 4 
gives an experimental test of this relation. For 
normal liquids, the experimental points fall well 
along a straight line: for associated liquids, more 
energy is required. As we have seen, for such 
liquids the heat of forming holes is no longer the 
determining factor, but rather the necessity of 
the flowing molecules breaking hydrogen bonds. 
An examination of the entropy of activation 
(Fig. 5) shows that these same associated liquids 
show a high positive A.S*; in terms of our model, 
after the flowing molecule has broken its hydro- 
gen bond it can rotate freely, and thus gains 
entropy. The entropy change tends to compen- 
sate the heat change,’ so that when AF* is 
plotted against AE,,, an excellent straight line 
is obtained (Fig. 6). The absolute value of the 
fluidity can be predicted, knowing only the 


4 A. E. Stearn and H. Eyring, J. Chem. Phys. 5, 113 
(1937). 


JOURNAL OF APPLIED PHYSICS 











molal volume and the heat of vaporization, from 
the equation : 


7) 


o=— exp (—AE yap/2.45RT). 


(22) 
Nh 


V. THe Mixture LAW FoR FLUIDITY 


The viscosity of a mixture of liquids is not 
related to the viscosities of the pure components 
by any simple additive relation such as 


n=Nimt+Non. (23a) 


Several equations have been tested in the search 
for an additive function for fluidity, among them 
being 


= Nigit Nodo, (23b) 
¢'= Nidii+No¢?!, (23c) 

o" = Nidi"+ Nod", (23d) 
log = N, log 1+ N2 log $2, (23e) 


the weighting being done according to weight 
fraction, volume fraction and mole fraction. There 
is also a wide variety of equations containing 
one or more adjustable constants. For example, 
an equation of the type of Eq. (23d) but with an 
adjustable exponent has been used by petroleum 
engineers to estimate the viscosity of mixtures 
of lubricating oils; for high-viscosity paraffin- 
base plus low-viscosity naphthene base, the 
exponent is 1/30; for the opposite case, the 
exponent is 1/3.1; for two oils of the same base, 
the exponent is 1/6.5.'6 Such equations are 
convenient for interpolation purposes, but that 
they fit the experimental data is more a matter 
of arithmetical inevitability than of merit. Of the 
equations not containing an adjustable constant, 
Eq. (23e) fits the experimental data rather 
better than the others. 

If the flow process were strictly determined by 
the properties of one molecule flowing, Eq. (23b) 
would be expected to hold. If it were determined 
by the properties of two molecules flowing past 
each other, Eq. (23c) would be expected to hold. 
However, if the presence of holes is necessary 





% Cf. J. Kendall and K. P. Monroe, J. Am. Chem. Soc. 
43, 115 (1921). 


16 G. Wilson, Nat. Petroleum News, May 22, 1929. 
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Squares: N2, CO, A, Oo, C2H«, HCl, HBr, HI, H2S. 


Circles: Typical organic liquids. 


Triangles: Water, alcohols, organic acids, glycerol. 
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for flow, the ability of the solvent to provide 
holes is the determining factor. Thus the average 
thermodynamic properties of the entire liquid 
will be involved when any individual molecule 
flows. 


As a simple approximation, Eq. (20) may be 
used for mixtures if for V is inserted the actual 
value of the average molal volume, and for AF* 
is inserted the weighted arithmetic mean of the 
values for the pure components: 

V 
o= 


Nh 


exp [—(NiAF+N24F2*)/RT]. (24) 


When V, and V2 are not too different, Eq. (24) 
reduces to Kendall’s equation, Eq. (23c). 

As a test of Eq. (24), AF* has been plotted 
against mole-fraction for a number of pairs of 
liquids. Three types of curves are obtained, of 
which typical examples are shown in Fig. 7: 

(a) Closely similar liquids, e.g., benzene and 
anisole, give a straight line. 

(b) Liquids which definitely form a compound, 
e.g., chloroform and ether, give a convex curve. 

(c) Liquids which are slightly dissimilar give 
a slightly concave curve; liquids which are 
markedly dissimilar, e.g., benzene and alcohol, 
give a markedly concave curve. 


In order to give a quantitative explanation of 





3.05 





2.55 











Fic. 7. 


676 





this effect, we recall the model of flow. When 
liquids are mixed, there is a free energy change 
(in addition to the change of ideal mixing) due 
to new bonds being made. When bonds are 
broken to make a hole, an additional amount of 
free energy is required which is a fraction of 
the free energy of this “‘quasi-reaction.’’ The 
deviations of AF* from a linear law (calories, at 
50 mole percent) have been plotted against the 
deviations from Raoult’s law (calories, at 50 
mole percent) in Fig. 8, for systems for which 
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1. Chloroform—ether 12. Acetone—ether 
2. Acetone—aniline 13. Benzene—aniline 
3. Chloroform—acetone 14. Benzene—cyclohexane 
4. Acetone—nitrobenzene 15. Ethyl alcohol—ether 
5. Ethyl alcohol—methy] alcohol 16. Benzene—cyclohexanol 
6. Benzene—toluene 17. Benzene—ethyl acetate 

Benzene—anisole 18. Carbon disulfide—acetone 

7. Benzene—ether 19. Ethyl alcohol—carbon tetra- 
8. Benzene—carbon tetrachloride chloride 
9. Benzene—nitrobenzene 20. Benzene—phenol 
10. Benzene—carbon disulfide 21. Benzene—ethyl alcohol 
11. Carbon disulfide—ether 22. Benzene—acetic acid 


partial-pressure data are available or can be 
estimated. It is seen that the experimental points 
do show a definite trend, and moreover that they 
tend to fall along a line drawn through the zero 
point with a slope of 1/2.45, as was found in 
Fig. 4. In view of the fact that these experi- 
mental points represent small deviations in large 
experimental quantities, the correlation is quite 
satisfying. The mixture law of Eq. (24) is to be 
modified to read 
Dd 


o= ——~ a0 |- (wnari+mars 


Nh 
AF® 
-) / 8th (25) 
2.45 


where AF® is the excess free energy of mixing.!’ 
~ 1G, Scatchard, J. Am. Chem. Soc. 53, 3186 (1931); ibid. 
56, 995 (1934); ibid. 57, 1805 (1935); Trans. Faraday Soc. 
33, 160 (1937). 
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The application of this mixture law is illus- 
trated in Fig. 9 where the experimental values for 
the fluidity of the system benzene-phenol are 
plotted, together with the curves calculated 
according to several different mixture laws. 


VI. Mixture LAW For DIFFUSION 


The absolute rate theory treatment of the self- 
diffusion of liquids has given the expression! 


kT 
D=)*?— exp (—AF*/RT) (26) 
1 
or, upon combining with Eq. (20), 
Dn=d kT, AoA3. (27) 


If we picture a solute, partitioned at equi- 
librium between two different solvents, it is clear 
that there is no net diffusion across the interface. 
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Hence it is the activity, and not the concentra- 
tion, which is the driving force in diffusion. 
The appropriate variants of Eqs. (26) and (27), 
which will be found useful in the treatment of 
diffusion in imperfect solutions, may be derived 
as follows :'8 

Consider a diffusing molecule passing over the 
potential barrier of Fig. 10, with the concentra- 





INITIAL ACTIVATED FINAL 
STATE STATE STATE 
l d | 

—— 
Concentration: et 
dx 
Activisy dy dy 
Coefficient: n nrg. NtAG, 
Fic. 10. 


tions and activity coefficients as illustrated. The 
general rate expression will be used in the form 


kT a. Yn 
k’ =— exp (—AF*/RT)—. (28) 
h 73 
Then the forward and backward rates of trans- 
port are 








kT 2 
adie exp (—AF*/RT) 
Y1 
— a 
addy, 
fue 
vi dx 
dni\ kT o 
rate, = (n+0- ‘), —exp (—AF*/RT) 
dx h 
q A dy1] 
mi it+t—— 
4 we dx J 
x—————_-_ (30) 
ad dy) 
vii 1+-— — 
be ¥Y1 dx J 


18 A. E. Stearn, E. M. Irish and H. Eyring, J. Phys. 
Chem. 44, 981 (1940). 
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Chloroform Ether 


Fic. 11. 


and the net rate is, upon simplifying, 


dn, kT o 
rat€net = ——? : -exp (—AF*/RT) 


dx 
d\n 7Y1 
x| 1+ | (31) 


Nn my, 








The diffusion constant D is defined by the 
relation 


dn, 
rate net = —D— (32) 
dx 
and upon combining (31) and (32), 
kT ‘ din 7 
D=)*— exp (— AF? Rn|i+ =| (33) 
h din n, 


If V; and V2 are not too different,* Eq. (33) 


* Strictly speaking, when the units of concentration are 
changed to mole-fraction, the small correction factor 
[1—N,(1—V,/V2)] must be introduced; however, since 
this is a linear function of mole-fraction, it will simply be 
absorbed into the corresponding values of 2’. 
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may be put in the form 
kT . ° dina, 
D=v*— exp (—AF* Rr) “| (34) 
h din N, 
Or the equivalent form 
NRT TF dina, 
Dn= -| — “| (35) 
AcA3 d In Ni 


As a test of Eq. (35) Dn and 


dina, 
of 
din N,; 


have been plotted against mole-fraction for the 
systems chloroform-acetone and chloroform-ether 
in Figs. 11 and 12.'* The resulting curves provide 
excellent confirmation of the equation. 


2.0 





1S 














Chloroform Acetone 


Fic. 12. 


Even when the vapor-pressure data is inade- 
quate for a quantitative check, Eq. (35) gives a 
good qualitative picture of the variation of 
diffusion with the composition. An example in 
the case of the system water-n propyl alcohol is 
shown in Fig. 13. 

Equation (35) was derived on the assumption 
that a single molecule is the diffusing species; 
if diffusion were a bimolecular process, the 
diffusion equation would be given by Eq. (20) of 


17H. Lemonde, Ann. de physique 9, 539 (1938). 
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reference 18. Thus the experimental evidence is 
that the flow process involves a single molecule 
passing over a barrier into a hole. 


VII. FLow ORIENTATION 


Absolute values of the diffusion coefficient 
predicted from Eq. (27) are consistently higher, 
by a factor of 2 to 5, than the values observed. 
This can be explained if the ratio A2A3/A1, is 
taken as greater than the average diameter, i.e., 
if the flowing molecule so orients itself that \, is 
its smallest dimension, and Xo; its largest area. 
Pure liquids in motion should then exhibit 
optical anisotropy, just as larger colloidal par- 
ticles do. This effect in pure liquids was dis- 
covered by Maxwell in 1873. 

We shall assume the following simple model 
for the Maxwell effect: If a molecule takes a 
particular orientation so that its area exposed to 
shear is \2’A;3’, its index of refraction is m’; the 
average molecule, with area 2°; exposed to 
shear, has an index of refraction n°. Then the 
two kinds of molecules will have different 
energies, and the fraction of oriented molecules 
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TABLE IV. 

SUBSTANCE M 1010 
o-dichlorbenzene 1.85 
p-xylene 1.84 
mesitylene 1.31 
m-xylene 1.29 
chlorbenzene 1.22 
o-xylene 1.21 
toluene 1.04 
phenyl ethanol 0.67 
benzene 0.64 
heptanol-1 0.41 
carbon tetrachloride 0.06 
cyclohexane 0.03 


will be given by the Boltzmann factor 


n’ 
—=exp [ —f(A2/A3’ —A2°A3°)x/RT | 
n°? 


=1—f(do/As’—2Az")x/RT, (36) 


which must be averaged over all x from zero to 
the top of the potential barrier. The quantity 
measured in the experiments, An/n, is given by 


— = ———————— = |] — — = f(do'd3' = A2°A3°) / 


+ 


r/2 
f x exp (—E/RT)dx 
0 











kT: = — (37) 
i) exp (— E/RT)dx 
0 
Maxwell’s experimental formula is 
An 
—=fM, (38) 
n 


where M is constant for a given substance. As an 
example, we may calculate M for benzene at 
300°K: do’A3’ is about 16X10-'* cm?, YA 2A°s 
11X10-'* cm?; the integration, carried out nu- 
merically for a 3-kcal. sinusoidal barrier, gives for 
the quantity in brackets (A/2)(0.145). The calcu- 
lated Maxwell constant is 0.52 10-". For com- 
parison, Table IV lists experimental values of the 
Maxwell constant for several liquids.” 

It is seen that Eq. (37) not only gives a func- 
tional dependence agreeing with experiment, but 
also predicts correct values for the constant. 


Physik 112, 407 (1939). 
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Correlation Between Elastic Moduli and Viscosity of Liquids and Plastics 


ANDREW GEMANT 
The Detroit Edison Company, Detroit, Michigan 


(Received May 21, 1941) 


The purpose of the paper is to indicate a method by means of which elastic constants (bulk 
and Young's moduli) and viscosity can be computed, one from the other, under conditions in 
which the direct determination of one of the data required is difficult or not feasible, and the 
other of the two quantities in question is known. 


1. INTRODUCTION 


HIS paper is concerned with a functional 

correlation between the elasticity moduli, 
including compressibility, of a liquid or plastic 
material and its constant of viscosity. 

The chief purpose of this presentation is a 
practical one. The correlation mentioned appears 
to be valid also for plastic materials or for a series 
of materials ranging from the liquid to the plastic 
state. It is often required to know viscosity 
values of plastics where direct measurements are 
difficult or even impossible, yet the measurement 
of an elastic modulus is easy. On the other hand, 
it is often required to know the velocity of 
propagation of acoustical waves (i.e., compressi- 
bility) of liquids for which viscosity measure- 
ments can be carried out much more easily and 
quickly. These requirements can be met by 
means of the above-mentioned correlation in the 
following way. 

The equation representing this correlation 
contains generally two parameters, their values 
depending upon the material in question and the 
temperature. Applied to a given substance, the 
equation connects compressibility and viscosity 
over the whole temperature or pressure range of 
the material. It further appears that for a given 
temperature the value of these parameters can be 
considered as constant, at least for practical 
purposes, even for a whole series of chemically 
related materials, by using a suitable average 
value of the parameter. If, now, for such a 
material or series of materials the two constants 
of the equation are known from a few available 
data, then the unknown data can be estimated 
from the value of their counterpart in a range 
where direct measurements are impractical. 

It should be added that the experimental 
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result of two papers,'? to be discussed below, 
have led the author to the conclusion that a 
correlation of the above type must exist. By 
means of quantitative data contained in one of 
the papers it was possible for the author to set up 
an equation on an empirical basis, connecting 
compressibility and viscosity. As a further step, 
theoretical equations referring to compressibility 
on the one hand and viscosity on the other were 
consulted, with the result that an approximately 
correct theoretical equation, relating the two 
data and having exactly the form of the above- 
mentioned empirically obtained formula, could 
be set up. The value of the parameters occurring 
in the equation could now be estimated theo- 
retically, and the agreement in the order of 
magnitude between computed and experimentally 
obtained values was quite satisfactory. Thus the 
significance of the empirical relation was con- 
siderably strengthened. 

In the following pages a brief theoretical 
deduction of the relation in question will be given 
and, in a number of cases for which numerical 
data are available, the relation will then be 
verified. The possibilities of practical applications 
previously referred to will be discussed along 
with the numerical cases. 


2. THEORY 


In seeking a relation between bulk modulus 
(reciprocal of compressibility «) and viscosity 7, 
the internal energy up» of 1 g of the material, as 
determined by intermolecular forces, will be 


1R. Matteson and C. J. Vogt, J. App. Phys. 11, 658 
(1940). 

2J. M. Davies and W. F. Busse, ‘‘Dielectric properties 
of plasticized polyvinyl chloride,”” Am. Chem. Soc. Meet- 
ing, Cincinnati, April, 1940; also Davies, Miller and 
Busse, J. Am. Chem. Soc. 63, 361 (1941). 
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considered. It is possible to establish a relation 
between x and up on the one hand, 7 and u“, on the 
other. Eliminating uo, the required relation x(n) 
is obtained. Again it should be stressed here that 
the following deduction is only approximate. 
However, it is felt that a more rigorous derivation 
will be worthwhile only if the suggested relation 
should really prove to be of practical value in the 
future. 

The relation x(uo) is based upon the existence 
of attractive and repulsive forces, as treated by 
M. Born, F. London, J. E. Lennard-Jones and 
others, and used for the explanation of the 
physical characteristics of liquids. The specific 
energy Uo is here considered as a sum of two 
terms, the repulsive and the attractive, both 
being certain inverse powers of the molecular 
distance. Let us denote this inverse power as 3n 
for the attractive and 3m for the repulsive forces. 
The value of m is around 3 to 4, that of m 1 to 2. 


It is then possible to show® the following relation 
to hold: 


k= —1/mnpto (1) 


with p=density of the material. This is the 
required relation «x(uo); it is an approximation, 
strictly valid only at absolute zero temperature. 

The next relation required is n(uo). There are 
several theories of viscosity,‘ but all agree in an 
equation of the form 


n=No exp (— lLMuo/RT). (2) 


Here mo is the product of certain molecular 
constants, / is a numerical factor, M=molecular 
weight, R=gas constant, 7=absolute tempera- 
ture; and the work against external pressure is 
neglected. By dividing by the density p, the same 
equation can be written down for the kinetic 
viscosity p: 
w=po exp (—1Muo/RT). (3) 
By eliminating u» from Eqs. (1) and (2), and 
writing 
A=0.43M1/mnRT (4) 
and 
B= —logio Mo, (5) 


3R. H. Fowler, Statistical Mechanics (Macmillan, New 
York, 1936), p. 140. 

‘J. M. Burgers, Second Report on Viscosity and Plasticity 
(Nordemann, New York, 1938), Chapter I. 
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one obtains 
A/x=p(logio u+B), (6) 


the desired equation correlating « and yu, and 
containing two parameters, A and B. 

It should be added that B itself is a function of 
the temperature and the density. The nature of 
this function depends upon the deduction of the 
theoretical equation (2) of which there are 
several in the literature. We do not wish to follow 
up this point here; it should be mentioned, how- 
ever, that for practical purposes it suffices for B 
to be considered either as a constant or as a 
quantity simply proportional to p, in which case 
one has, instead of (6): 


A/x=p(logio u+Cp). (7) 


Equation (6) and, in one particular case, (7) will 
be used in the next sections. 

The possible practical applications of Eq. (6), 
as pointed out above, are twofold. First, it can be 
used for one given substance, in order to correlate 
x and y for varying temperature or pressure. The 
energy uo changes with both of these variables in 
an unknown manner; there is, however, no need 
to know this, since “» does not explicitly enter the 
equation. Second, the formula can be used for a 
series of related materials. Theoretically, both A 
and B do change from one member to another of 
such a series, but for practical purposes an 
average value assumed for A and B appears to 
give sufficiently accurate values in calculating 7 
from «x, or vice versa, as the case may be. 


3. APPLICATIONS 
(a) Hydrocarbons 


As a first example, heptane in a temperature 
range from 10° to 60°C should be considered. 
Equation (6) will be used. The constant A will be 
defined as referring to room temperature (75 in 
deg. K), and thus the equation for other tempera- 
tures T has to be written: 


To 1 
A-—-—=p(logio n+ B). (8) 

7 a 
Used in this sense, the equation takes the form of 
the well-known viscosity-temperature relation 
with the important difference that the factor of 
1/T now varies with the temperature, as does 1/x. 
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In other words, the slope of the log u(1/7) plot is 
not a constant, but decreases with increasing 
temperature, since the bulk modulus also de- 
creases with increasing temperature. Such a 
deviation from the straight line is frequently 
observed. Whereas, as said above, Eq. (6), and 
thus (8), are by no means theoretically accurate, 
they are so for practical purposes, and so the 
deviation from the straight line of the log u(1/T) 
plot is also sufficiently explained on this basis. 
We shall come back to this point when discussing 
oils. 

The data for heptane are shown in Table I. 
Column 1 gives the temperature in degrees 
centigrade, column 2 the density and column 3 
the viscosity. The next two columns give values 
of the bulk modulus, column 4 those determined 
by means of measurements? of sound propagation 
velocity, and column 5 those calculated from 
Eq. (8) with the following two constants: 


A=0.031X10-°; B=2.63. 


As can be seen, the agreement is satisfactory. 

Next a series of various hydrocarbons is con- 
sidered. Since the constants in this case are just 
average values, the formula is limited to related 
substances, as will be seen from Table II in which 
column 1 gives the name of the material, column 
2 the density, column 3 the viscosity at room 
temperature, columns 4 and 5 the bulk modulus 
at room temperature as measured? and as calcu- 
lated from Eq. (6) with 


A=0.091X10-°; B=3.43. 


Whereas the agreement is good for the saturated 
members, it is less satisfactory for tetralin, con- 
taining double bonds. This indicates that if (6) is 
to be used for a series of materials, A and B are 
constants only as long as related substances are 
considered. 

Let us now check the empirical values of A and 
B in the light of Eqs. (4) and (5), deduced 
theoretically. It will be seen that the order of 
magnitude agrees very well. In order to get a 
plausible value for the hydrocarbon series, let us 
take hexane as representative of the series. For 
evaluating B, one of the existing viscosity 


5S. Parthasarathy, Proc. Ind. Acad. Sci. 2A, 497 (1935). 
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TEMPERA- n, CENTI- 1/«x, ABS. UNITS 

TURE, °C r) POISES MEASURED CALCULATED 
10 0.692 0.465 9.92 x 10° 9.9 10° 
20 0.684 0.416 9.14 9.1 
30 0.675 0.375 8.37 8.3 
40 0.667 0.341 7.65 7.8 
50 0.658 0.311 6.94 7.0 
60 0.649 0.284 6.34 6.4 


TABLE II. Computations on various hydrocarbons (20-25°C) 
by means of Eq. (6). 


n, CENTI- 1/x, ABS. UNITS 

SUBSTANCE p POISES MEASURED CALCULATED 

Pentane 0.626 0.235 6.90 x 10° 6.9 10° 
Hexane 0.660 0.304 8.15 7.9 
Heptane 0.684 0.403 9.07 9.0 
Octane 0.704 0.520 10.0 10.0 
Cyclohexane 0.779 0.93 12.4 12.8 
Tetralin 0.971 2.00 20.7 18.4 





theories has to be consulted, and accordingly the 
most recent one, that by Eyring® should be 
selected, according to which 


no=1.1X10°M'T!/VIE 


(M=molecular weight, V=molar_ volume, 
E=heat of evaporation per mole). For hexane 
M=86, V=130 cc, E=6220 cal. Thus pyo=0.49 
xX 10-* and B=3.3, which agrees with the value 
assumed to satisfy the experimental data. In (4) 
we have the factor / which, according to Eyring, 
is about }. For n the value 2 is most likely ; that 
for m is uncertain, but is quite probably 3. Thus 
A=0.06X10-°, again in satisfactory agreement 
with the assumed value. 


(b) Petroleum oils 


One of the papers, mentioned in the Intro- 
duction, is by Matteson and Vogt! on the velocity 
of acoustical waves in petroleum fractions. They 
found that the velocity was a function of the 
density and viscosity of the oil. This finding gave 
the author the idea that a functional correlation 
between x and 7 must exist. 

The experimental data are best reproduced by 
Eq. (7) in this case with 


A=0.45X10-*; C=10.5. 


¢R. H. Ewell, J. App. Phys. 9, 252 (1938). 
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The result over a number of oils of a fairly 
extensive series is shown in Fig. 1. 

The figure shows the sound propagation 
velocity ({[1/px ]*) as a function of the kinematic 
viscosity, the density of the oils being the 
parameter of the set. The agreement between 
individual curves seems only moderate, par- 
ticularly for the highest viscosities. From what 
was said before, this cannot be expected to be 
otherwise. Besides, one has to consider that it is 
not one curve, but a whole set of curves that is 
here reproduced by an equation containing two 
constants only. 

The usefulness of the equation in this case is 
evident, and follows from the paper of Matteson 
and Vogt. For any given petroleum fraction, to be 
used in Diesel engines, only density and viscosity 
have to be determined (both easy, routine 
operations), and the velocity of compressional 
waves, the knowledge of which is important for 
their operation, can be computed with fair 
accuracy. 

For a given oil the relation can be used to 
estimate viscosity or compressibility for any 
temperature or pressure if one of the two data is 
known. 
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Let us check the validity of the relation in this 
respect. As to the variation with temperature, 
Eq. (8) has to be used as pointed out before. In 
other words the log n(1/T) plot should not have a 
constant slope, as generally assumed in a first 
approximation, but should vary as 1/x—in a 
second approximation. Looking at a _ typical 
transformer oil between —20° and +80°C,’ the 
slope at —20° to that at +80° is in a ratio of 


7 Andrew Gemant, Liquid Dielectrics (Wiley, New York, 
1933), p. 26. 
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about 2 : 1. As to x, data from the International 
Critical Tables* for kerosene between +20° and 
+80°C, and extrapolated to —20°C, can be 
taken. The extreme values for 1/x are: 16.4 and 
9.6 10°, and their ratio is 1.7. Thus a large part 
of the variation of the log n(1/T) slope can be 
accounted for by the corresponding change of x, 
indicating the usefulness of the relation. 

As to the pressure effect, it is known that both 
the bulk modulus and the viscosity increase with 
pressure. Moreover, the increase of the viscosity 
is exponential, that of the modulus linear, in 
agreement with Eq. (6). But we can go a step 
further, and see how far the numerical agreement 
goes. The viscosities of a mineral oil for 1100 and 
2200 atmos. (and 40°C) are given by R. B. Dow 
and co-workers? as 1.5 and 9.0 poises. Writing 


1/x=1/kxo+bP, (9) 


where }=a constant and P pressure (in atmos.), 
and substituting into Eq. (7), b can be computed 
to 1.310°*. In the Critical Tables,* b for kerosene 
and some lubricating oils is given as around 
4X 10°, which agrees, as to its order of magnitude, 
with the value computed from viscosities and 
using Eq. (6) or (7). This choice makes no 
difference here since p is practically independent 
of the pressure. 


(c) Polyvinyl chlorides 


The second paper referred to in the Intro- 
duction is one by Davies and Busse? on the 
dielectric properties of plasticized polyvinyl 
chlorides. The authors quoted found pronounced 
power factor maxima with varying temperature, 
the curves shifting towards lower temperatures 
with increasing tricresyl phosphate content. Such 
curves show that the opposing forces of dipolar 
rotation are of a viscous nature, in contradistinc- 
tion to other solids in which it is the solid friction 
that has to be overcome by the dipoles. 

Now the maximum of power factor curves 
occur when the time constant of the dipoles 
coincides with the reciprocal of the a.c. frequency. 
Since the time constant according to Debye is 
proportional to the viscosity, it follows that the 
differently plasticized polymers at the respective 








8 Int. Crit. Tab. 2, 146 (1927). 
®R. M. Dibert, R. B. Dow and C. E. Fink, J. App. 
Phys. 10, 113 (1939). 
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temperatures of maximum power factor must 
possess the same viscosity. The viscosity of such 
polymers is not easy to measure and, instead, 
hardness figures were obtained by Davies and 
Busse by means of the Shore Durometer, with 
the surprising result that these figures at the 
temperatures in question were very nearly the 
same. 

The Durometer figure is chiefly a somewhat 
involved function of Young’s modulus of the 
material ; in other words, the elastic modulus was 
the same at the temperatures of power factor 
maxima. This result is understandable only if it is 
assumed that there is a definite relationship be- 
tween the viscosity of the polymers and their 
elastic moduli. This conclusion is in complete 
agreement with the theory developed in these 
pages. The data presented by Davies and Busse 
are not of a nature to allow a check of the 
equation, as was possible for the simple hydro- 
carbons and the petroleum fractions. 


(d) Glasses 


Glasses are among the objects believed to be 
particularly suited for being applied in con- 
nection with Eq. (6). The reason is that the 
elastic modulus of glass (Young’s modulus E in 
this case) is easily measured at lower tempera- 
tures; its viscosity, however, only at higher 
temperatures. The missing data, often quite 
important, can then be computed if the constants 
are known. 

A check of Eq. (6) in this case is not easy and 
can be only approximate. The reason for this is 
that there are no corresponding data for both 7 
and E obtained on one particular glass by any one 
author. The data had to be compiled from work 
of different authors on similar, but certainly not 
identical, glasses. First we replace 1/x by E, 
using the relation 


1 pE 
alii daniiimmnciietapii (10) 
k (1+p)(1—2p) 


where p= Poisson’s ratio, in the average 0.25 for 
glass. Thus, 1/x=#/2.5. Further, p can be 
replaced by its average value 2.5, and its varia- 
tion neglected. One then has for Eq. (8): 


110AE/T=2.5 log n+(2.5B—1). (11) 





TABLE III. Glass at different temperatures. 


TEMPERA- LOG10 7 

TURE, °C E iN ABS, UNITS EXPERIMENTAL COMPUTED 
400 0.57 x 10" 17 17.6 
450 0.55 15.8 15.4 
500 0.53 13.7 13.4 
550 0.51 11.4 11.6 
600 0.48 9.8 9.8 


TABLE IV. Glass of different sodium content at 600°C. 


PERCENT Ex 

Na2O LOGi0 9 EXPERIMENTAL COMPUTED 
25 8.7 0.57 x 10" 0.57 X 10" 
23 9.3 0.61 0.60 
21 10.0 0.63 0.63 
19 10.4 0.66 0.65 
17 10.7 0.67 0.66 
16 11.0 0.68 0.68 
15 11.4 0.69 0.69 
14 11.8 0.71 0.71 


All data for the computations given in the 
following tables were taken from the monograph 
on glass by G. W. Morey.” 

In Table III data are given for glass of about 
the following composition: SiOQ,: 68 percent, 
Na2O+K,0: 20 percent, Ca0+ MgO: 8 percent, 
Al.O3: 4 percent. Young’s modulus E was 
measured by Badger and Silverman from room 
temperature to 450°C, and the data up to 600°C 
are extrapolated values. Viscosity data were 
obtained by Scott, Irvine and Turner. The 
constants in Eq. (11) for this case are 


A=0.60X10-, B=5.2. 


The agreement between calculated and de- 
termined 7 values is satisfactory, and the order of 
magnitude of the constants is the same as in the 
previous cases, their absolute values necessarily 
being different. 

Table IV summarizes data for glasses of 
different composition, all at the same tempera- 
ture of 600°C. The SiO, content was around 74 
percent, and the Na,O content is shown in the 
table. Viscosity values are taken from measure- 
ments after English. Corresponding data of E are 
available only for room temperature (after Clarke 
and Turner), and these are shown in column 3. 
From the previous series at different tempera- 


10 G. W. Morey, The Properties of Glass (Reinhold, New 
York, 1938). 
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tures it could be concluded that the ratio Ego0/E20 
is about 0.74. Consequently the E data obtained 
from Eq. (11) were divided by 0.74 in order to 
bring them to room temperature level. The 
constants of the equation used for this case were 


A=0.58X10; B=4.1 


and the calculated E data are shown in the last 
column, showing good agreement with the 
experimental values. Hence, the usefulness of the 
formula for glass too is shown. 


SUMMARY 


By means of a simple theory an equation, 
establishing a connection between elastic moduli 
of liquids and plastics and their viscosity, is 
derived. The equation, which contains two con- 
stants only, can be used either for a given 
substance in a range of varying temperature and 
pressure, or for a series of chemically related 
materials. The usefulness of the equation is 
shown on simple hydrocarbons, petroleum frac- 
tions, plasticized polyvinyl chlorides, and glasses. 





The Solution of A.C. Circuit Problems 


Louts A. PIPEs 
Graduate School of Engineering, Harvard University, Cambridge, Massachusetts 


(Received May 31, 1941) 


Starting with the set of linear equations that determine the complex currents of a general 
n-mesh circuit, a method is presented for the numerical solution of the equations to determine 


the various mesh currents of the system. 


The method presented has computational advantages over the conventional methods usually 
used in that only real numbers are involved in the various numerical operations and only the 
evaluation of second-order determinants is involved. The operations required are those of matrix 
multiplication and are most conveniently performed with a calculating machine. The paper is 
concluded by the solution of a general Y—Y circuit with a neutral return as an example of the 


method presented. 


INTRODUCTION 


ANY articles have appeared in the litera- 

ture during the past few vears on the 
application of matrix algebra to electrical prob- 
lems. In all these papers, the concept of the 
inverse of a square matrix is involved. Since in 
the solution of alternating current problems, the 
elements of the various matrices are complex 
numbers, it is seen that the amount of compu- 
tational labor involved in the inversion of these 
matrices is enormous. 

For example, consider a general four-mesh 
network. To determine the four-mesh currents 
from the impressed complex voltages and the 
complex impedances of the system we must 
invert a square matrix of the fourth order whose 
elements are complex numbers. This involves 
the evaluation of the determinant of the matrix 
itself which, in this case, is one of the fourth 
order and also the evaluation of sixteen determi- 
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nants of the third order if the original matrix is 
not symmetric or ten third-order determinants 
if it is symmetric. Since in general all these 
determinants have complex elements, it is seen 
that the labor involved is very great. 

In this paper a method is presented that 
attacks the problem in an entirely different 
manner. By the use of this method, the compu- 
tational labor is greatly reduced since only real 
quantities are involved and the procedure is one 
of matrix multiplication which is so easily per- 
formed by means of a calculating machine. The 
elements of matrix algebra presented in (1) will 
be assumed in this discussion. 


1. THE GENERAL EQUATIONS 


Let us consider a general m-mesh circuit in the 
steady state. Let the various complex mesh 
currents have the values (11, Js, ---J,) and let 
the complex voltages impressed on the contours 
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of the meshes be (£;, Eo, ---E,). The voltages 
and currents are related by the equations: 


ZilitZilet+ ees +Zinl,=Ei 


Jrdeitoveih cies 10 
“LnwdlitZnelet+ sale Nn +Zanln= E, 


where the complex numbers Z,, are the various 
self-impedances of the meshes and Z,, where 
r#s are complex numbers that are the mutual 
impedances of the system. In static networks, 
the condition: 


) a) (2) 
holds. 
The set of linear Eqs. (1), may be written in 
the compact matrix form: 


[Z}(1) =(E), (3) 


where [Z] is a square matrix of the mth order 
whose elements are the complex numbers Z,,, 
(J) is a column matrix whose elements are the 
various complex mesh currents and (E£) is a 
column matrix whose elements are the complex 
mesh voltages. 

In the usual problem of circuit analysis, the 
mesh voltages and the complex impedances are 
given and it is required to determine the complex 
mesh currents. The fundamental problem is 
therefore to solve the set of Eqs. (1) for the 
complex currents, [,. 


2. THE ADMITTANCE MATRIX 


If the circuit under consideration is a static 
one, then the condition (2) holds and the matrix 
[Z] is a symmetric matrix so that we have: 


(ZY =(Z], (4) 


where [Z]’ is the transposed matrix of the 
matrix [Z] obtained by interchanging the rows 
and the columns of [Z]. In the general case, 
the matrix [Z_] may be written in the form: 


(2]=(RJ+jLX], (5) 


where j=1/—1, [R] is a square matrix of the 
nth order whose elements are the various self- 
and mutual resistance coefficients of the circuit, 
and may be termed the resistance matrix. [X ] 
is a square matrix of the mth order whose 
elements are the various self- and mutual re- 
actance coefficients of the system. 
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Let us now define a matrix, [Y], by the 
equation: 


CYJ=(Z), (6) 


where [Z}! is the inverse of the impedance 
matrix [Z]. This matrix, [Y], is called the 
admittance matrix of the system. The elements 
of [ Y ] are given explicitly by: 


Y,,= Cofactor of Z,,/|Z\, (7) 


where |Z| is the determinant of the matrix, 
[Z]. Since [Y] is the inverse of [Z], it has the 


property that: 


[ZILYJ=LVILZJ=U, (8) 


where U is the unit matrix of the mth order. 
From (7), it is seen that if [Z] is a symmetric 
matrix, then [Y] is also a symmetric matrix 
and satisfies the relation: 


LYY=LY], (8a) 


where [Y]' is the transposed matrix of [Y]. 
In view of (6), if we premultiply both sides of 
(3) by LY], we obtain explicitly: 


(1)=(Z}"(E) =LY](£). (9) 


In (9) we thus have an explicit expression of 
the various mesh currents in terms of [Z }-! and 
the mesh voltages. In most articles on the 
application of matrix algebra to circuit theory, 
the Eq. (9) is given as the solution of the problem 
under consideration. However, when one comes 
to the actual numerical solution of a problem of 
this type, it is seen that the amount of computa- 
tional effort required is enormous. 

We now turn to the development of a pro- 
cedure which materially reduces this computa- 
tional labor and allows us to work with real 
numbers entirely. 

Let us write: 


CLYJ=(G)+j[BJ=(2Z)", (10) 


that is, we separate the complex matrix [Y ] 
into its real and imaginary parts. The square 
matrix of the mth order, [G] will be called the 
conductance matrix and the square matrix of 
the nth order, [ B ] will be called the susceptance 
matrix. From (7), it is seen that if [Z] is a 
symmetric matrix, then both [G] and [B] are 
symmetric matrices. We may now perform the 
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multiplication : 


(ZILVYJ=(Z][ZJ"'=U 
=(LRJ+jLX])((G]+jB)) 
= (LRILG]—LX IB) 
+j(CRILBI+LX IG). 


Equating the real and imaginary coefficients of 
this expression, we obtain: 


U=(RI](G]—LX IB), 
LOJ=(X ][GJ+LRILB], 


where U is the unit matrix of the mth order and 
[0] is the zero matrix of the mth order. These 
two equations may be solved for the matrices 
(G] and [B] and we obtain: 


(11) 


(12) 
(13) 


(GJ=(LRJI+LX IRI OD), (14) 
(BJ=—-LR}'LX ILC] (15) 

" (B]=—-((XJ+(RIEX}P LR), (16) 
LG]=—-([X]}"LRILB]. (17) 


Either the first or the second set of equations 
may be used to compute [G] and [B]. Which 
set is used depends on the complexity of the 
matrices [X ] and [R]. 

The computation of [G] and [B] by either 
of the two sets of equations is greatly facilitated 
by the use of a powerful iterative method of 
determining the inverse of a non-singular matrix 
whose elements are real numbers. This method 
will be described in a later section. 


3. DETERMINATION OF THE MESH CURRENTS 


To determine the mesh currents having 
obtained the conductance and_ susceptance 
matrices, let us decompose the voltage and 
current matrices into their real and imaginary 
parts as follows: 


(E) = (c) +j(d) 
(I) = (a) +3(0). 


(18) 
and 


(19) 
We then have: 


(1) =LYJLEJ=((6)+j LB) LO +310]. 


Equating real and imaginary coefficients, we 
obtain: 


(20) 


(a)=[G](c)-[BI@), 
(6) =(B](c)+[G](d). 


(21) 
(22) 
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The elements of the column matrices (a) and 
(b) give the real and imaginary parts respectively 
of the various mesh currents. 


4. PARTITIONED MATRICES 


Before proceeding to the discussion of a 
powerful iterative method for the inversion of 
real non-singular matrices that enables the 
operations indicated by the set of Eqs. (14) and 
(15) or (16) and (17) to be performed, it is 
necessary to digress briefly to discuss the subject 
of partitioned matrices since the method to be 
presented makes use of this concept. 


(a) The order of a matrix 


A matrix is said to be of the (m, n) order if it 
has m rows and nm columns. The first quantity, 
m, refers to the rows, and the second quantity, 
n, refers to the columns. 


(b) Multiplication of matrices of different orders 


Two matrices of different orders can be 
multiplied together only if a certain condition is 
satisfied. If the number of columns of a matrix, 
[ B] is equal to the number of rows of a matrix, 
[A], then these two matrices are described as 
conformable. They can then be multiplied in the 
order [BA]. If [B] is a matrix of the order 
(q,m) and [A] is a matrix of the order (n, p), 
then the product [B][A] is a matrix of the 
order (q, p). This rule may be expressed sym- 
bolically in the form: 


(q, n) X(n, p) =(q, p). (23) 
For example, 
4 2 -1 2 2 3} 
3-7 1 -—8||-3 0 
2 4 -—3 1 1 
> & 
ae 9 
=| 4 6}. (24) 
[-8 -8 





(c) Partitioning a matrix 


It is possible to divide the array of a given 
matrix into smaller arrays by horizontal and 
vertical lines. The matrix is then said to be 


687 








partitioned into submatrices. As an example, 


consider 
1 9;3 
[a]=|2 8)4 (25) 
6 2,7 


each of the sub- 
matrices of a partitioned matrix by a single 
symbol and then the original matrix becomes a 
matrix of matrices. 


It is possible to represent 


Partitioned matrices may 
be multiplied as if the submatrices were ordinary 
matrix elements provided they are conformable. 
An examination of this criterion leads to the 
rule that in general for the multiplication of two 
partitioned matrices it is necessary and sufficient 
that for every partitioning line between columns 
of the matrix on the left there shall be a parti- 
tioning line between the corresponding rows of 
the matrix on the right and no further horizontal 
partitioning. For example: 


ea; | Gee | Aas bo boo bos 


ayy adie fi bie bis 


31 Q32 | 33 bs bso bss 


Qi, Qy2 | Aig] foun bie bis 
=|d21 G22 | de3|}be1 dee bes}. (26) 


31 32 | A33'bbs, D323 


THE INVERSION OF MATRICES 
»9F SUBMATRICES 


Many methods for the inversion of non- 
singular matrices may be devised. (See reference 
2.) Of these methods, the one making use of 
submatrices appears to lend itself most readily 
to numerical computation. This method will 
now be described. 

Let us suppose that we have given a non- 
‘singular matrix, [A |, of the mth order, and we 
wish to determine its inverse, [B]. We then 
must have: 


BY THE USE 


[B][A]= Un, (27) 


where U,, is the unit matrix of the mth order. 
Now let the matrices [A ] and [B ] be partitioned 
in the following manner: 
A 12(7, S) 
(28) 


A 40(s, S)_ 


Axl, r) 


A — 
[4] r) 
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“Bilr, r) 


B\= 
LB] ~ r) 


B21, S) 
| (29) 


Boo(s, S) 


The orders of the submatrices are indicated 
by the notation of Section 4, and we have 


r+s=m. Performing the indicated operation 
(27), we obtain the relations: 
ByAutBypAn= U,, a 
ByAi2t+BypA»=(0], 31) 
BA 1+Bx2A 21 =([0], 32) 
BA 12+ B2A 22= Us, (33) 


where U, and U, are unit matrices of the rth 
and sth orders, respectively. These four equa- 
tions may now be solved to give the four sub- 
matrices of [B]. If we use the notation: 


Masi Av. (34) 
V=AnAn, (35) 
W =Ao— VAi2=A2x2—AaM, (36) 


we then have: 


Bu=An+MW V, (37) 
By = —MW-, (38) 
Bu = —W-"V, (39) 
By =W-". (40) 


These equations determine the required matrix 
. . —l —1 
[B] provided that the inverses A,;,; and W 
P , ak ve , 
exist. When A, is found, the quantities M, V, 
W can be calculated and [B] determined. The 
rumerical work may be arranged in the follow- 
ing convenient tabular form : 














4 21 Los 
M=A Ae A il A 12 
w- V=AnAy W=An— VA 
then 
11 o* MW-'V | 


|\—-MW- 
— |. (41) 


A [4 = 
LA Loy Ww 


If we choose r=m—1, s=1, then A»; is a single 
row, Ai is a single column, 
element, and W is a scalar. 


Ax is a single 
This suggests the 
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foll 
nu: 
un 


29) 


ted 


ave 


10n 


30) 
31) 
32) 
33) 
rth 


la- 
ib- 


9) 





following procedure for dealing with a given 
numerical matrix: Suppose that the matrix 
under consideration is: 


"An Ai. +++ Ain 
[A]=]An Aso ++: Arnl. (42) 
Ant An. — ; A * 


We may then set up the following sequence of 
matrices, S;, Se, S3, etc. which progressively 
include more and more of the principal diagonal 
elements of [A ] formed in the following manner: 


S:=[(Au], (43) 
A A 2 

s.=| * (44) 
Ax Ao» 
Ay Aw Ais i So | Ais 

S3= = A 2 Ao _— | Ags ’ (45) 
A31 A 32 A33 As, Azo | Ass 


[Ai, Ai. Ais Ais 
Ax, Azzy Ass Ads 
As: Az. Asz3 Ass 
[Aa Age Aas Aas 


S4 








| Ais) 
- S Ax . (46) 
: | Asa \ 
ha tes di ted 


Since Ss is a second-order matrix, its inverse is 
easily computed and is given by: 


im 1 A2 —A1s 
5° — re | | (47) 
(Ay Ax2—AwAn)L—-An An 


; , 1 

Now using this value of S: , we can now apply 

to S; the scheme of computation described above 
7 = 

and thus obtain S; . In the same manner, we use 

‘ 1 . 

S; and apply the process to Sy. In this way we 

‘ ‘ a 

finally obtained S, or [A }". 

In exceptional cases it is possible that one 
may encounter a singular matrix, say the S; 
matrix during the iterative process. This is 
indicated by obtaining W=0 for that particular 
matrix. In such cases, the simplest procedure is 
to transfer the ith row of [A ] to the last place 
and to continue in the normal manner. This 
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Fic. 1. General Y—Y circuit. 


rearrangement of the rows is equivalent to a 
premultiplication of [A] by a matrix N, where 
N is the unit matrix of the mth order with its 
ith row transferred to the last place. ‘ 

The final reciprocal matrix will then have to 
be post-multiplied by N to give [A]. The 
above iterative process for: the computation of 
the reciprocal of a non-singular matrix is very 
easily carried out with a computing machine 
since the principal process required is that of 
matrix multiplication. 


6. NUMERICAL EXAMPLE 


To illustrate the above theory, let us apply 
it to the calculation of the performance of the 
general Y—Y circuit of Fig. 1. For generality, 
we shall consider mutual impedances between 
the three sending-end impedances, the four line 
impedances, and the three receiving-end im- 
pedances. The sending-end potentials assumed 
are quite unbalanced. The problem is to compute 
the three line currents J,, J, and J, and the 
current in the neutral. 


Assumed numerical values 


The following numerical values are assumed: 


Zaa=10+ 74 La=Zva=Jj2 Zy=2-j5 


Zw=3-j7 = Zue=Zu=j3 Z2=3+j8 
Zee=4tj9  LZea=Lac=jl Zss=1+j6 
Zi2=Zn=j4 Zr=2+j6 Zae=Ze=j3 
Zis=Zs=j2  Zu=84+j3 Zm=Zyr= je 
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Z32=Z23= 3 Z1=6+j2 


Znn=9+j3 


Zn=Ziu=j4 
Zen=Zns= Jl 
Zu =Zu=j3 
E,.=110+ 70 Zur=Zin=j2 
E,= —58— 749 
E.= —5+ 115. 
By writing the Kirchhoff second law and 
taking the potential drops around the three 


independent meshes including the neutral return 
of the system, we obtain: 


“Gir Aye Ay3) TL E, 
G21 Ge d23||J,)}=| A} =[(Z](J). (48) 
Q31 G32 334% I. E.- 


The a;; quantities have the values: 

€11 = Laat Zrr— 2Zar+Zant+Ziu= 23+ j4, 
A22=Zowt+Zee— 2Znst+Zunt+Z22= 234+ J5, 

433 = Leet Zit— 2Znit+Zant+Z33 = 20+ j16, 
12=Zavt+Zre—Zin— Znst+ Zan tZi2=an=9+j7, 
O13 = Lact Zri—Zin— ZnttZantZi3=431= 9+ j6, 
423 = Zoet+Zet— Zen — Zntt+ Zant Z23=432=9+79. 


Separating the matrix [Z] into its real and 
imaginary parts, we obtain: 


23 9 9 4 7 6 
crRJ=| 9 23 9|, (x]J=]7 5. 91. 
9 9 20 6 9 16 


In this example, the resistance matrix is particu- 
larly simple and it is more expedient to calculate 
the conductance and susceptance matrices by 
using the Eqs. (14) and (15) rather than (16) 
and (17). Before computing the inverse of [R] 
we may first divide out the factor 9 and write: 


2.555 1 1 
[R]=9| 1 2.555 1 |=9[A] 
1 1 2.222 
and hence : 
2.555 1,17} 
[R}'=3] 1 2.555; 1 |=s{[A}. 
eo 
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| To find the inverse of the above square matrix, 
we may now use the procedure developed in 
Section 5. First we compute the inverse of S, 
as follows: 
5 2.555 cal 
“11.000 2.555)’ 
0.46228 Spree 


| — 0.18089 0.46228 


II 


_— 
Se 





From this, we now set up a table as in Section 5 
as follows: 


1.00000 1.00000! 2.222 | 
M 0.28139 0.46228 —0.18089| 1,000 
0.28139 | —0.18089 0.46228! 1.000 


0.60261 | 0.28139 0.28139! 1.659. 
w V Ww 

0.047714 0.047714 

MW-1V= a 
0.047714 0.047714 


Substituting into (42) of Section 5, we obtain: 


0.50999 —0.13318 —0.16956 
[A}'=| —0.13318 0.50999 —0.16956 
—0.16956 —0.16956 0.60261 
and hence: 
([R}'=3[4T' 
~ 0.05666 —0.01479 —0.01884 
=} —0.01479 0.05666 —0.01884]. 
—0.01884 —0.01884 0.06695 
We now compute: 
0.01007 0.22442 0.19446 
| [R}'LXJ=] 0.15311 0.01021 0.37647 
—0.09459 0.11976 0.78860 
and 
| CXICRT'LX] 
2.77798 2.94273 5.19156 
=| 2.94273 4.51105 7.03407}, 
5.19156 7.03407 13.12790 


| CRI+CXILRTI LX] 
| 25.778 11.942 | 14.192 

=|11.942 27.511 | 16.034| =[A]. 
| 14.192 16.034 ; 33.128. 
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We now must compute the inverse of this new 
matrix [A]. We choose S: and compute its 
inverse as follows: 





25.778 11.942 
S2=! 11.042 meal 
_, [ 0.04855 —0.02107 
52 =] _ 9.02107 canal 


Again we compile the table: 
| 14.192 16.034 | 33.128 


m | 0.35119 | 0.04855 —0.02107| 14.192 | 
| 0.43036 | —0.02107 0.04549| 16.034 


0.04707 | 0.35119  0.43036| 21.244 














| 





w V W 

0.00580 0.00711 

MW" v-| | 
0.00711 0.00871 


Hence substituting into (41) of Section 5, we 
obtain: 


(LRI+LX IRI TX =[4T? 


0.05435 —0.01396 —0.01653 
=| —0.01396 0.05420 —0.02025| =[G]. 
—0.01653 —0.02025 0.04707 


We therefore have obtained the conductance 
matrix [G], and note that it is symmetric. 
From (15), we now obtain the susceptance 
matrix, [B] as follows: 


[B]=—-(RI"LX IG] 


0.00028 —0.01007 0.00772 
=} —0.01007 0.00500 —0.00172}. 
0.00772 —0.00172 —0.03712 


We note that both the matrices [G] and [B] are 
symmetric matrices as they should be. This is 
a check on the arithmetical work. To obtain 
the mesh currents, we use Eqs. (21) and (22) 
and obtain: 
5.488 
(a) =[G ](c) —[B](d) =| —4.135], 


4.176. 
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1.342 
(6) =[B](c) —LG](d) = | — 6.373 
7.539 
since in this case, 
110 —. 
(c)=|—58], (d)=| —49}. 
—5 115 


The mesh currents are, accordingly: 


T,=5.488+ 71.342, I,= —4.135—j6.373, 
T.=4.176+j7.539. 


The current in the neutral is: 
T,= —(Uatlo+J-) = — 5.529 — 72.508. 


The problem is thus completely solved. 


7. CONCLUSION 


From the above example, it is seen that the 
amount of computational labor involved in the 
numerical solution of a.c. problems by the 
method presented in this paper is substantially 
less than that involved in the conventional 
procedure. The principal type of numerical 
operation involved is that of matrix multiplica- 
tion. This operation is most simply performed 
by the use of a calculating machine. The pro- 
cedure avoids operations with complex numbers 
entirely, and the work is extremely systematized. 
The saving in labor is more apparent in circuits 
having a large number of meshes. If the network 
under consideration has only resistance, or only 
reactance coefficients, then it is only necessary 
to find the inverse of [R] or of [X ] respectively. 
This inversion may be carried out conveniently 
by the method of submatrices. 
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Surface Studies with the Electron Microscope 
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P to the present time the electron microscope : 


has achieved its greatest successes in the 
observation of objects thin enough to transmit 
readily electrons incident with 50 to 100 kv 
velocity—objects, like bacteria, viruses, smokes, 
thin films, which are materially less than 1 
micron in thickness. The direct methods of 
observing surfaces of compact bodies have not, so 
far, led to comparable image quality at high 
magnifications. Among these methods may be 
mentioned: ° 

(1) the observation of surfaces with the aid of 
thermionic electrons, photo-electrons, or second- 
ary electrons emitted from the surface with very 
small velocities,'~* and 

(2) surface observation with reflected electrons; 
here the surface may either be normal to the 
viewing direction,® the incident electrons striking 
it at a glancing angle, or it may be nearly parallel 
to both the direction of observation and the 
direction of illumination.® In the first case the 
wide velocity distribution of the reflected elec- 
trons makes high resolution impossible, while in 
the second case extreme foreshortening and 
unequal resolution in two mutually perpendicular 
directions are serious drawbacks. 

Another method designed to accomplish the 
same purpose consists of scanning a small area of 
the surface with a very fine electron probe and 
reproducing the surface on the screen of a 
cathode-ray tube or some other type of recorder 
by the method of wired television.7~* It is 
‘impossible to gauge, at the present time, the full 
range of application of this last method. Results 


1'E. Briiche and H. Johannson, Naturwiss. 20, 353 
(1932). 

2M. Knoll, F. G. Houtermans, and W. Schulze, Zeits. f. 
Physik 78, 340 (1932). 

* E. Briiche, Zeits. f. Physik 86, 448 (1933). 

*V. K. Zworykin, J. Frank. Inst. 215, 535 (1933). 

®°E. Ruska and H. O. Miiller, Zeits. f. Physik 116, 366 
(1940). 

* B. v. Borries, Zeits. f. Physik 116, 370 (1940). 

7M. Knoll, Zeits. f. tech. Physik 16, 489 (1935). 

*M. v. Ardenne, Zeits. f. Physik 109, 553 (1938). 

°F. R. Moulton, Sci. Monthly 49, 189 (August, 1939). 
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published so far, while very promising, have not 
materially exceeded the capabilities of the light 
microscope. However, active work which may 
lead to greatly improved performance is now in 
progress. 

In view of the limitations of the direct methods 
of surface observation, it was decided to investi- 
gate the possibilities of adapting the standard 
transmission-type electron microscope’! to sur- 
face studies by preparing suitable replicas of the 
original surface, these replicas being thin enough 
to readily transmit 50-kev electrons. 

Such a procedure had already been suggested 
by H. Mahl, who prepared thin oxide films on 
aluminum surfaces"-" and floated them off after 
etching away the metal underneath with the aid 
of a saturated mercuric chloride solution. 

Mahl also mentioned the possibility of 
generalizing the method by using collodion films 
in place of oxide films, and prepared some replicas 
of aluminum surfaces using the aforementioned 
method for removing the film from the metal." In 
addition to destroying the specimen, this method 
has the drawback of leaving large quantities of 
impurities on the replica unless great care is taken 
in washing it. It is, furthermore, restricted to 
materials which can be dissolved without 
affecting the replica film, i.e., primarily metals. 

A large number of different methods of pre- 
paring replicas have been investigated by the 
authors. One method in particular has given 
excellent results, appears most general in its 
application, and does not destroy the specimen. 
Two alternative methods, which will be described 
later, share the advantage of not damaging the 
surface if properly applied. 

A thick layer of metal, e.g., silver, is evapo- 
rated in a vacuum chamber onto the surface to be 
studied. This metal film is stripped mechanically 

10. Marton, Phys. Rev. 58, 57 (1940). 

1 J. Hillier and A. W. Vance, Proc. I. R. E. (in print). 
2H. Mahl, Zeits. f. tech. Physik 21, 17 (1940). 

18H. Mahl, Jahrbuch d. AEG-Forschung 17, 67 (1940). 


4H. Mahl, Metallwirtschaft 19, 488 (1940). 
1H. Mahl, Metallwirtschaft 19, 1082 (1940). 


JOURNAL OF APPLIED PHYSICS 





er 


re- 


en 
its 
on. 
ed 
he 


0- 
be 


lly 


[CS 








Fic. 1. Bainite (0.78 percent carbon steel). Positive Fic. 2. Fine Pearlite (0.78 percent carbon steel). Positive 
replica obtained by silver-collodion process. Magnification: replica obtained by silver collodion process. Magnification : 
20,000. 20,000. 




















Fic. 3. Glazed porcelain surface. Positive replica obtained Fic. 4. Quartz cleavage surface. Positive replica obtained 
by silver-collodion process. Magnification: 20,000. by silver-collodion process. Magnification: 20,000. 
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from the surface, yielding a negative metal 
replica of the surface. If too thin for the stripping 
process, the metal film may be reinforced by 
electroplating. Next, a dilute solution of the final 
replica material, e.g., a 1 percent solution of 
collodion in amyl acetate, is flowed over the 
negative replica surface and is permitted to dry. 
Then the metal replica with the adhering thin 
film is immersed for several hours in a suitable 
solvent for the metal—2 to 3 normal nitric acid 
in the case of silver. At the end of this period the 
metal is completely dissolved and, after washing 
in distilled water, the residual positive replica 
film may be placed on the fine-mesh screen (e.g., 
250-mesh copper or stainless steel) which acts as 
the object support in the electron microscope. 
The choice of metal to be used as transfer 
medium is dictated primarily by the adhesion 
properties of the evaporated metal with respect 
to the surface to be studied, and by the ease of 
evaporating and dissolving it. Silver has been 





Fic. 5. 
‘periodic error. Positive replica obtained by gelatin-col- 
lodion process. Magnification: 8000. 


15,000 line/inch diffraction grating with large 


found very satisfactory in the case of steels, 
brass, aluminum, various alloys, and non-metallic 
substances. On the other hand, for a silver surface 
some other material, such as aluminum, would be 
required. 

As a first example of the application of the 
method just described, Fig. 1 shows an electron 
micrograph of Bainite (suitably heat-treated 0.78 
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percent carbon steel).* The same steel, with 
different heat treatment, has the fine pearlite 
structure seen in Fig. 2. Figures 3 and 4 show 
examples of the application of the method to non- 
metallic surfaces, i.e., a glazed porcelain surface 
and a cleavage surface of natural quartz, re- 
spectively. The heavy black particles visible in 
Figs. 3 and 4 are silver residues, which aid 
materially in the focusing of the image in the 
electron microscope. 

Two other alternative methods which, though 
more limited in application, have given some 
interesting results, merit brief description. The 
first of these employs gelatin{ in place of evapo- 
rated metal as the transfer medium. A concen- 
trated aqueous solution of gelatin is spread on the 
surface of the specimen and is permitted to dry 
thoroughly. The gelatin will then either separate 
itself from the surface or can be split off from the 
surface with the aid of a razor blade. The 
solution of final replica material is flowed on the 
negative gelatin replica so obtained, dried, and 
the gelatin dissolved in a warm dilute solution of 
ammonium hydroxide. After washing with dis- 
tilled water, the residual positive thin-film replica 
is transferred to the object-supporting mesh. 
This method is limited to objects which are not 
affected by the water in the gelatin solution, and 
thus cannot be used for the study of ordinary 
steels, for example. Furthermore, it was found 
that, with quartz and glass surfaces, the adhesion 
of the dried gelatin is often so strong that it 
cannot be removed without destroying part of 
the surface, limiting the usefulness of the method 
also in this case. As an example of the application 
of the gelatin-collodion method, Fig. 5 shows an 
electron micrograph obtained from a metallic 
copy of a 15,000 lines per inch grating with a very 
large periodic error. 

The other alternative method to be discussed 
leads to a negative rather than a positive replica 
of the surface. It consists of direct mechanical 
stripping of the thin film to be observed from the 
surface of the specimen. The dilute solution of 





* The metallographic specimens used in this study were 
prepared and supplied to us by Dr. R. F. Mehl, Director 
of the Metals Research Laboratory of the Carnegie Insti- 
tute of Technology in Pittsburgh. Dr. Mehl’s cooperation 
in this respect was an essential factor in obtaining the 
results here reported. 

+t Method suggested by Dr. L. Marton. 
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Fic. 6. Bainite (0.78 percent carbon steel). Negative 
collodion replica stripped from specimen. Magnification: 
20,000. 


replica material is flowed directly on the specimen 
surface and dried. A more concentrated solution 
is applied near one edge of the surface, forming a 
much thicker film at this point. This thick film, 
after drying, is raised with the aid of a razor 
blade. Using this raised portion as a handle, the 
film, or as much of it as will adhere to the thicker 
portion, is stripped off and placed on a clean 
water surface, from which it is finally transferred 
to the object-supporting mesh. Both collodion 
and vinylite lacquer* have proved to be suitable 
replica materials for this procedure. While ex- 
tremely simple, it is relatively uncertain, and 
yields far poorer resolution than the first method 
discussed. Figures 6 and 7, obtained from the 
same Bainite and Pearlite specimens as Figs. 1 
and 2, demonstrate its limitations. Finally, the 
possibility of greatly distorting the extremely 
thin film during stripping cannot be disregarded. 


* Kindly furnished us by Dr. Mehl. 
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Fic. 7. Fine Pearlite (0.78 percent carbon steel). Nega- 
tive vinylite lacquer replica stripped from specimen. Mag- 
nification: 20,000. 


In conclusion, it is worth while to point out one 
basic advantage of the replica methods as com- 
pared to the direct methods of surface observa- 
tion with the electron microscope. The density at 
any point of the picture of a positive replica is a 
direct measure of the height of that point with 
respect to a reference plane parallel to the sur- 
face. Ridges in the surface show as dark portions, 
depressions as bright portions of the image. 
Furthermore, if the exposure and the develop- 
ment procedure are standardized, image density 
can be correlated once and for all with replica 
thickness, so that the actual contour of any cross 
section of the surface imaged can be derived 
directly from the corresponding microphotometer 
record of the negative. The microscope is 
calibrated in advance by taking pictures either of 
superposed collodion films of known thickness or 
of evaporated films of known mass density. 
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Silver-Magnesium Alloy as a Secondary Electron Emitting Material 


V. K. Zworykin, J. E. Ruepy, anp E. W. PiKe* 
Research Laboratories, RCA Manufacturing Co., Inc., Camden, N. J. 
(Received May 1, 1941) 


HE wide use of secondary electron emission 

electrodes as a means of obtaining gain in 
thermionic vacuum tubes has been retarded by 
the lack of a sufficiently rugged emissive ma- 
terial. There are many such electrodes that can 
be built into tubes which will function satis- 
factorily for a few minutes, or a few hours, but 
the requirements are for thousands of hours of 
useful life. The failure with time of these initially 
satisfactory emitters is usually due to a combi- 
nation of the effects of electron bombardment, 
residual gas, and temperature rise. The operating 
conditions are different for nearly every type of 
tube, but, roughly, the minimum requirements 
for a satisfactory secondary emission electrode 
are that it shall be able to operate at 300°C with a 
primary current density of 5 milliamperes per 
square centimeter in a reasonable vacuum for 
1000 hours, with not more than 20 percent change 
in gain. 

It had been known for a number of years that 
certain compounds of magnesium, when suitably 
combined with a conducting material, gave high 
secondary emission, but not generally 
satisfactory in operating stability. On the basis of 
this information, some metal alloys containing 
magnesium were obtained and tested. The results 
were so promising that an investigation was 
undertaken to determine the most suitable alloy 
and the best technique for handling it. It was 
found that while alloys of a small percentage of 
magnesium with copper, silver, gold, and alumi- 
num all had essentially the same secondary 
emission properties, the addition to some of these 
of a few percent of calcium, titanium, or beryllium 
did not produce any marked difference. Further 
work with alloys of different constituents was 
then discontinued in order to investigate the 
properties of the most promising of this group. 
The amount of magnesium in the alloy was found 
not critical between less than 1 percent and 15 
percent, except in the matter of ductility. Various 
methods of making it, by both air and vacuum 


were 


* At present with Andrews and Petrillo, Inc. 
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melting, gave satisfactory results. The simplest 
procedure was to melt a known weight of silver in 
a steel crucible in an electrically heated oven, 
quickly thrust a piece of magnesium beneath the 
surface of the melt by means of a steel wire, stir 
a little, and then allow the mixture to cool. The 
resulting ingots were sawed and rolled to a 
convenient thickness for the cutting out of 
electrodes. 

After an electrode was formed, a final cleaning 
of the active surface was carried out in any one of 
a variety of ways, such as by washing with soap 
and water, alcohol, or acetone, or by scraping, 
sanding, buffing, or electrolytic etching. The 
electrode was then mounted in the desired 
structure and sealed into a glass tube, care being 
taken to avoid excessive heating of the alloy in 
the process. 

When such a tube was simply pumped without 
baking or other special treatment, a gain of about 
1.5 was observed. (Unless otherwise stated, all 
values of gain given are for the case of 200-volt 
primary electrons.) If, however, the tube was 
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Fic. 1. Gain vs. primary electron energy for activated alloy. 


baked as in a normal outgassing routine, the gain 
was usually between 4 and 6. Induction heating 
of the alloy in the presence of oxygen would 
further raise the gain by several points. The gain 
as a function of bombarding voltage for a typical 
activation of this type is shown in Fig. 1. 
Accompanying the increase in gain, the alloy 
surface changed from a bright silver to a golden 
or yellow color. Examined by means of electron 


JOURNAL OF APPLIED PHYSICS 





St 


1g 


It 


iS 





diffraction, unactivated alloy showed only the 
pattern of silver, and, after activation, only the 
pattern of magnesium oxide. Since various prepa- 
rations of magnesium oxide were known to have 
comparable gains, it was assumed that activation 
of the alloy was simply due to the formation of 
some of this compound, and that its superior life 
properties lay in the better binding of the oxide to 
the metal base. 

The successful activation of alloy surfaces was 
a difficult process. Aside from those treatments 
which obviously affected the surface, and which 
might have been expected to influence the gain, 
there were some obscure causes of variability. A 
considerable amount of work was done in in- 
vestigating these, and the conclusion was finally 
reached that minute traces of some gases or 
vapors given off from the inside of the tube 
during the initial bake were the principal factors 
in the uncontrolled variability. The alloy itself 
was found to be as potent a contamination carrier 
as any other part. The speed with which the 
contaminant was removed from the neighborhood 
of the active surfaces was also important, and for 
this reason different activation techniques were 
required for tubes of different geometries. 

Satisfactory life characteristics are a matter of 
prime importance; hence life tests were under- 
taken as a routine procedure. A schematic cross 
section of the elements of the tube used for this 
purpose is shown in Fig. 2, together with its 
associated circuit. W is a ring of 0.005-inch 
tungsten wire used as cathode, mounted inside 
and connected to the nickel cylinder K in such a 
position that the alloy target A is not exposed in 
a direct line to any part of W. C is the nickel 
collector disk. The primary current J, is adjust- 
able by means of the rheostat R. The gain is the 
ratio I,/I,. The primary electrons strike the 
target over an area slightly greater than 3 inch in 
diameter. Since the distribution over the bom- 
barded area is not uniform, the observed life 
curves (gain vs. time) are composites, covering a 
range of current densities. 

During early study of the alloy, the activation 
aim was always to obtain maximum initial gain, 
after which the tubes were sealed off, and placed 
on life test at an average primary current density 
of about 5 milliamperes per square centimeter. 
During the runs some of these became definitely 
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gassy and exhibited rapid drops in gain, while the 
best ones showed a rather rapid decrease during 
the first fifty hours, and then a gradual flattening 
out to a slow linear decrease. The typical life 
curve of a high initial gain surface is shown in 
Fig. 3. Some surfaces which had been in operation 
for a few hundred hours and had decreased in 
gain were rejuvenated by the admission of a little 
oxygen, but after about thirty more hours of 
operation the gains were back on their original 


Ar 
2v ac. FR! 
Ww. © a 
Ip \~ 


ie) C 





“ 























-_—7 
Ov ie 1% 
3 nag 
+400v 











Fic. 2. Schematic section of life test tube elements and 
associated circuit. 
curves as if there had been no reactivation at all. 

The simplest procedure used for obtaining the 
type of activation described above was as follows: 
The tube was baked until outgassed. After it had 
cooled, the alloy and other metal parts were 
outgassed with high frequency. Oxygen was then 
introduced and the alloy heated by induction, 
after which the oxygen was pumped out, the 
alloy and metal parts again heated with high 
frequency, the Ba-Al getter exploded, and the 
tube sealed off. 

The hope was held for some time that it would 
be possible, by proper pumping and activation, to 
produce surfaces that would maintain their high 
initial gains, but after building more than fifty 
tubes and always finding essentially the same life 
characteristics, it was decided that the decay 
effect was something fundamental. 

Attention was then turned to the activation 
processes which gave lower initial gains, to see if 
some of them would yield surfaces of greater 
stability. The typical life curve of a surface 
activated by a modification of the procedure 
described above is shown in Fig. 4. The gain is 
here quite good, and the stability excellent after a 
preliminary aging. Another property of the alloy, 
which makes it valuable for use in high gain 
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Fic. 3. Alloy activated for maximum initial gain. 


multipliers, is its very low dark current,! as 
compared with the commonly used types of 
emitters employing alkali their 
formation. 

Although evidence upon which to base an 
explanation of the secondary emission behavior of 
the alloy is not conclusive, one reasonable 
hypothesis is that the stable gain of about 4 is 
due to the silver-magnesium oxide itself, while 
the higher gains are due to oxygen adsorbed or 
absorbed in the silver-magnesium surface. The 
initial rapid decrease in gain with time of 
operation is due primarily to the surface losing its 
absorbed oxygen. Thereafter the behavior is a 
function of ion bombardment, loss of oxygen by 


~ 1Z. Bay, Rev. Sci. Inst. 12, 127 (1941). 
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Fic. 4. Alloy activated for maximum stable gain. 


the surface, and rejuvenation by oxygen from 
other parts of the tube. 

Another possible explanation is that the high 
gain is due to the Ag-MgO and that the decay 
under bombardment is caused by the decompo- 
sition of the MgO with the liberation of oxygen 
(which is taken up by the getter), and thus the 
gradual formation of a surface of silver mag- 
nesium. Curves of the type shown in Fig. 4 would 
be explained by the formation of a less de- 
composable silver-magnesium oxide structure. 

Although there are many points of both 
theoretical and practical interest that still need 
clarifying, the use of the alloy as a secondary 
electron emitter in certain types of thermionic 
tubes is now feasible. 





Glass-to-Metal Seals. II 


A. W. Hutt, E. E. BurGer, AND L. NAvIAS 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received June 25, 1941) 


This is a continuation of the work reported in 1934. The greater part of the paper is devoted 
to a discussion of the properties of iron-nickel-cobalt alloys, especially the composition Fe 54 
percent, Ni 31 percent, Co 15 percent—called fernico—including sensitiveness to impurities and 
polymorphism. Three new glass-metal combinations are described, the metals being 42 percent 
nickel-iron, 26 percent chromium-iron, and pure iron, respectively. The strength of the glass- 


metal bond is discussed briefly. 


1. INTRODUCTION 


N a previous article' we showed how the stress 

in wire seals can be calculated from the 
difference in expansion of wire and glass; and 
the calculated stresses in test seals were com- 
pared with photoelastically observed stresses 
for a number of common sealing materials. 
A new alloy of composition Fe 54 percent, Ni 28, 


1A. W. Hull and E. E. Burger, Physics 5, 384 (1934). 


698 


Co 18 was described, whose thermal expansion 
matched very closely, over the whole range from 
room temperature to the sealing temperatures, 
the expansion of a special glass, 705-AO, de- 
veloped by the Corning Glass Works for this 
purpose. It was shown that seals made with 
this alloy and 705-AO glass were practically free 
from strain. 

The present article gives the results of a 
continuation of these studies. It includes the 
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properties of iron-nickel-cobalt alloys, especially 
an alloy of composition Fe 54 percent, Ni 31, 
Co 15; which will be referred to as fernico. Two 
new glasses that seal to 42 percent nickel-iron 
and pure iron, respectively, are described. Some 
data on the strength of seals are given. 


2. FERNICO 


Iron-nickel-cobalt alloys have two properties 
which recommend them for sealing purposes: 
(1) Their thermal expansion in the neighborhood 
of room temperature is less than that of any 
other available metals or alloys, with the excep- 
tion of tungsten and molybdenum. This is im- 
portant for seals, because it allows the use of 
low coefficient glasses which are less likely to be 
broken by sudden changes of temperature than 
glasses of high expansion. (2) The thermal 
expansion of these alloys is non-linear like that 
of glasses, thus affording the possibility of a glass 
and metal with identical thermal expansion. The 
non-linear property of Fe-Ni-Co alloys appears 
to be associated with ferromagnetism; the ex- 
pansion being low in the temperature range in 
which the alloy is ferromagnetic, and increasing 
by a factor of three at the temperature (Curie 
point) at which ferromagnetism is lost. 

The expansion curves of the Fe-Ni-Co alloys 
which are of interest are shown in Fig. 1, com- 
pared to molybdenum, against which they were 
measured. It will be noted that the temperature 
of the transition from low to high expansion rises 
with increasing cobalt content, up to about 18 
percent. However, as this percentage is ap- 
proached, the alloy becomes less and less stable 
with respect to the y—a transformation (see 
next section). The composition 54 percent Fe, 
28 Ni, 18 Co represents the highest transition 
temperature for stable low coefficient alloys. 

The transition temperature of low expansion 
glasses, on the other hand, is in general high— 
higher than that of any of these alloys, or any 
known low expansion alloys; and the glasses 
become less stable, especially with respect to 
water solubility, the Jower their transition tem- 
perature. 

For strain-free seals the transitions of glass 
and metal should coincide, as well as their 
average expansions. This requirement could be 
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satisfied only by choosing a metal alloy having 
the highest possible transition temperature, and 
a glass with the lowest possible transition tem- 
perature. The combination of 705-AO glass and 
54 Fe, 28 Ni, 18 Co, given in our former paper, 
appeared to be the closest match obtainable. 

Since that paper was published, we have 
succeeded in matching the expansion and transi- 
tion of 705-AO glass with a still more stable 
alloy, viz., 54 Fe, 31 Ni, 15 Co. This alloy is a 
correct match only when free from all impurities, 
including manganese and silicon, which are con- 
sidered essential to the ductility of cast nickel 
alloys. We have found that pure Fe-Ni-Co 
alloys, without manganese or any other added 
element, made in hydrogen by the sintering 
process under proper conditions, are ductile and 
free from embrittlement under all the conditions 
of ordinary use. These conditions include heating 
and annealing in air, hydrogen, or gas; soldering ; 
brazing; and welding. The following Table I 
gives some of the physical properties of pure 
54 Fe, 31 Ni, 15 Co, made in this way. This 
alloy will be designated as fernico. 

Figure 2 shows the thermal contraction of 
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Fic. 1. Thermal expansion of iron-nickel-cobalt alloys, 
compared to molybdenum. 
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TABLE |. Physical properties of fernico: Fe 54, Ni 31, Co 15 


Tensile strength —72,000 to 80,000 Ib. per sq. in. 


Percent elongation —25 to 33 
Percent reduction in area —62 to 63 
Elastic modulus —18x 108 


Density 

Electrical resistivity 

Temperature coefficient of 
expansion —4.95 x 10~* from 25° 

Melting point —1460°C 

Hardness —75 ‘‘Rockwell B” 

Yield point —55,000 Ib. per sq. in 


—8.24 grams per cc 
—43.8 microhms per cm cube 


to 300°C 





fernico and of 705-AO glass. These curves were 
taken during cooling at the rate of 1 degree /min. 
from above the annealing temperature of the 
glass. Hence the close agreement of the curves 
from the annealing temperature down to room 
temperature means that glass and metal, if 
joined together at the annealing temperature, 
will contract at the same rate as they cool, and 
be free from strain both during cooling and after 
reaching room temperature. 

The confirmation of this prediction is shown 
by the insert in Fig. 2, which gives the photo- 
elastic stress pattern’ of a test seal made with 
these materials and cooled at the rate of 1 degree / 
min. In this pattern, stress in the glass is indi- 
cated quantitatively by bending of the inter- 
ference lines where they cross the samples. For 
the seal here shown, viz., 705-AO glass of 
diameter 7.5 mm on a fernico wire of 2.5 mm 
diameter, a deviation of one-tenth the distance 
between fringes would represent an axial stress 
of 0.20 kg/mm?. The lines in Fig. 2 are straight 
within one one-hundredth of a fringe, showing 
that the stress is less than 0.02 kg/mm‘. 


Dependence of stress on composition of alloy 


The chief criterion for freedom from stress in 
seals is close agreement of both the expansion 
coefficients and the transition temperatures of 
metal and glass. This agreement depends prin- 
cipally upon composition of glass and alloy. 
The sensitiveness of this agreement to composi- 
tion of the alloy, and hence the maximum 
allowable variations in composition, can be cal- 
culated easily, if expansion and transition tem- 
perature are expressed as functions of composi- 

2A. W. 


Hull and E. E. Burger, Rev 
(1936). 


7, Sci. Inst. 7, 98 
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tion.’ Instead of transition temperature, it is 
more useful to calculate the ‘‘sealing tempera- 
ture,”’ defined as the temperature of intersection 
with the expansion curve of a line 15 percent 
greater in slope than the average expansion 
coefficient between 0° and 300°C (see Fig. 2). 
A point so determined on glass curves represents 
quite accurately the sealing temperature of boro- 
silicate glasses used with these alloys; hence it 
also represents the sealing temperature of 
metal which matches the glass. 

Values of sealing temperature of iron-nickel- 
cobalt alloys as defined above were obtained 
from thermal expansion curves like those in 
Figs. 1 and 10, and are plotted in Fig. 3 as a 
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Fic. 2. Thermal contraction of fernico and 705-AO 
glass. The insert shows the stress pattern of a test seal 
between these materials. 


function of nickel and cobalt content. These 
data are represented satisfactorily over the com- 
plete range of low expansion alloys by a linear 
relation between sealing temperature and _ per- 
centage composition : 


T,=36.0 (Ni+0.875 Co—31.0) deg. C (1) 


when Ni, Co represent percent weights of Ni, Co. 

Similarly, the average expansion coefficient 
between 25°C and the sealing temperature is 
given by 


a,=0.805 (Ni+0.75 Co—34.9) X10-®. (2) 


These equations can be used to calculate the 


3H. Scott (Trans. A.I1.M.E., Inst. of Metals Div., 1930, 
p. 506) has given approximate equations for these quan- 
tities. 
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changes in length at sealing caused by variations 
57, and da, in sealing temperature and expansion 
coefficient, respectively, of the metal. It can be 
shown easily (see Fig. 4) that the differences in 
length at the sealing temperature of the glass 
between metal and glass (assumed equal in 
length at room temperature) due to variations 
67, and da, are: 


bly = — (tan a’+tan a,)57,, (3) 


51. = (T./cos* as) das, (4) 
where 7, is the sealing temperature, a, the 
average expansion coefficient between room tem- 
perature and sealing temperature, and a’ the 
coefficient (slope) at the sealing temperature. 

Since the angles are small, Eqs. (3) and (4) 
may be written: 

ily = = (a’ —a,)6T,, (3’) 

dla = T bas. (4’) 


Insertion of 57, and 6a, from (1) and (2) in 
Eqs. (3’) and (4’) gives: 
bly = 36a,(6 Ni+0.8756 Co), (5) 
5], =0.805 K 10-*®7 (6 Ni+0.756 Co). (6) 
For fernico, 
a’ =12.310-* 
a,=5.92«10-6 
T,.=472°C. 
Hence, 


ily = — 230 X10-*(6 Ni+0.8756 Co) 
dl. = 380 X 10-*(6 Ni+0.756 Co). 


—" 


Sealing Temperature 





——+ Percent Nickel 


Fic. 3. Sealing temperature of iron-nickel-cobalt alloys as 
function of nickel and cobalt content. 
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Fic. 4. Plot showing changes in the length of the metal 
at sealing caused by variations 57, and da, in sealing tem- 
perature and average expansion, respectively. Curve A is 
the thermal expansion curve of a correctly matched metal ; 
Curve B, a mismatched one. 
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Fic. 5. Thermal expansion cycle of iron-nickel-cobalt 
alloy with y—a@ transformation at 0°C, compared ta molyb- 
denum. Note the difference in expansion of the a and y 
phases. 


The total variation in length at sealing is 


6l = dlr + dla 
= 10~-*(1506 Ni+846 Co). (7) 


This difference 6/ in length at sealing is equal 
to the difference in contraction during cooling, 
and is a measure of the resulting stress. It was 
shown in our former paper,‘ that, for thick glass, 
a difference in contraction of 10-* cm/cm will 
cause a tensile stress of 0.55 kg/cm?. Hence the 
limiting stress of 100 kg/cm’, which was given 
in our former paper as the “breaking stress,” 
and which represents well our experience with 
the breakage of seals during the last five years, 
corresponds to a difference in contraction of 
100/0.55 X 10-* = 180 K 10-*cm/cm. Equation (7) 


4 See reference 1, Table I. e 
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shows that this difference in contraction, and 
hence a breaking stress, can be caused by a 
variation of 1.2 percent in Ni or 2.1 percent in 
Co. A seal stressed to the breaking point would 
not be satisfactory. If we define a satisfactory 
seal as one which has not more than one-fourth 
the breaking stress, then a total variation in Ni 
of +0.3 percent, or a total variation in Co of 
+0.5 percent (assuming no variation in Ni), 
could be tolerated. This is probably within the 
reach of good metallurgical practice. 

The situation is very different with regard to 
impurities. Only ‘‘desirable’’ impurities will be 
considered, viz., Mn, Si, and Al, the first being 
useful as a de-sulphurizer, the latter two as de- 


oxidizers. The equations for sealing temperature 
and average expansion at the sealing temperature 
of alloys containing these impurities are: 


T,=36.0(Ni+0.875 Co—0.5 Mn 


—1.4Si—1.0 Al—31.0), (1’) 
a, = 0.805(Ni+0.75 Co+0.5 Mn 
+0.6 Si+0.8 Al—34.9)x10-* (2°) 


and the changes in length caused by variations 
are: 
bly = — 230 XK 10-*(6 Ni+0.8756 Co 
| —0.55 Mn—1.46 Si—1.06 Al) 
380 X 10—-*(6 Ni+0.756 Co 
+0.55 Mn+0.66 Si+0.86 Al). 


bla 





FiG. 6. X-ray diffraction patterns of annealed fernico: (A) before and (B) after 4 hours in liquid air. Both show the face- 
centered lattice of the y phase, with no sign of transformation. 


(A) 








Fic. 7. Photoelastic 
stress patterns of fernico- 
705-AO seals: (A) before 
and (B) after 15 hours in 
liquid air. No strain has 
been produced by the low 
temperature treatment. 





Fic. 8. (A) X-ray diffraction of filings from annealed fernico rod. The body-centered lattice proves transformation 
to a phase caused by plastic deformation. (B) Pattern of same filings after 5 min. anneal in H2 at 900°C, showing restora- 


tion of y phase. 
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The total change, 6/= 6l7+6lq is: 


61 = 10-*(1506 Ni +845 Co+3056 Mn 
+5506 Si+5346 Al). (7’) 


Since breaking stress will result when 6/=180 
<10-*, the equation for breaking stress is 


180 = 1506 Ni+846 Co+3056 Mn 
+5506 Si+534 Al. (8) 


This shows that a breaking stress will result 
if the alloy contains 0.21 percent more or less of 
Mn and of [Si+Al ] than the prescribed amount ; 
while for a ‘‘satisfactory’”’ seal (<} breaking 
stress) the content of Mn and Si must each be 
accurate to 0.05 percent. 

The uncertainty in Mn and Si concentration 
in commercial alloys is generally much larger 
than this. 


Polymorphism of fernico 


The low thermal expansion of fernico is charac- 
teristic of its y (face-centered cubic) phase. 
This phase is metastable below about 400°C, the 
stable form being the a (body-centered) phase. 
The @ phase has a much higher thermal expan- 
sion than the y, as can be seen in Fig. 5, which 
shows the contraction and expansion, as the 
temperature is changed from +500°C to — 100°C 
and back again, of an alloy whose y—a trans- 
formation temperature is 0°C. The composition 
of this alloy, by analysis, is: 


Composition of Fe-Ni-Co Alloy No. 37518 
Fe—55.04 (by difference) percent 


Ni—27.86 
Co—16.59 
Si— 0.23 
Mn— 0.18 
C— .05 
P—<.02 
S— <.03. 


It is evident from Fig. 5 that the presence of 
even a small amount of the a phase would cause 
large stresses, and must be avoided. 

In the case of 54 percent Fe, 31 Ni, 15 Co 
fernico, the stability of the y phase is so great 
that the y—a transformation point, if it exists, 
is below the temperature of liquid air; i.e., the 
transformation to the body-centered form cannot 
be induced by temperature at any temperature 
down to and including that of liquid air. This is 
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Fic. 9. Thermal expansion of iron-nickel alloys. The 
expansion curves of the common commercial glasses are 
shown for comparison. 


1 


TT CONTRACTION OF Ni-Fe ALLOY 
+—}- __ AND GLASS No.1075 
_ | } + il - 
| 


Contraction in (cm per cm )xi07® 





300 400 500 
Temperature in Degrees C 


Fic. 10. Thermal expansion curves of 42 alloy and 1075 
glass, and photoelastic stress pattern of a seal between 
them. Absence of stress in the seal is shown by the straight- 
ness of the interference lines crossing it. 


shown in Fig. 6, in which (A) is the x-ray 
diffraction pattern of a thin strip of fernico 
annealed at 1100°C in hydrogen for one hour, 
and (B) that of the same strip after four hours at 
liquid-air temperature. Both show the perfect 
face-centered lattice of the y phase. 

The immunity to transformation at low tem- 
perature is also demonstrated by the strain 
patterns of a group of fernico-705-AO seals, 
shown in Fig. 7, (A) before, and (B) after 15 
hours in liquid air. Not only are the seals un- 
injured, but no strains can be detected as a 
result of the prolonged exposure to liquid air. 

While these tests prove that the y—a trans- 
formation of fernico (Fe 54, Ni 31, Co 15) cannot 
be induced by temperature, it may be produced 
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by plastic deformation, for example, by rolling, 
drawing, or filing. Fig. 8(A) shows the x-ray 
diffraction pattern of filings from a fernico rod 
which had previously been annealed. The filings 
show a mixture of the face-centered pattern and 
the body-centered cubic pattern of the a phase, 
caused by the plastic deformation of filing. 
After annealing 5 min. at 900°C, these filings 
gave the face-centered cubic pattern shown in 
Fig. 8(B), proving complete restoration of the 
vy phase. The y phase can be restored by a brief 
anneal, e.g., 5 minutes, at any temperature 
between 900° and 1200°C. But a longer anneal, 
viz., 4 hours at 900°C or 1 hour at 1100°C in 
hydrogen, is recommended in order to remove 
carbon, which causes bubbles in seals by reducing 
the oxides of the glass. 


3. NICKEL-IRON ALLOYS 


We have made a careful study of the nickel- 
iron system, in the hope of finding a combination 
of glass and nickel-iron alloy whose thermal 
expansions are “‘matched”’ in the sense defined in 
the last section. Figure 9, which gives the thermal 
expansions of the iron-nickel alloys in the low 
expansion range that is of interest, shows that 
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Fic. 11. Thermal expansions of 286 glass and 26 chrome- 
iron, and photoelastic stress pattern of a seal between them. 
Absence of stress is shown by the straightness of the dark 
lines crossing the seal. 
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the transition temperatures of these alloys are 
all much lower than those of commercial glasses. 
Hence a new glass had to be found to meet these 
conditions. This glass, whose General Electric 
designation is ‘1075 glass,” is a lead-borosilicate, 
of the following composition (see Table IT): 


TABLE II. Composition of 1075 glass. 


SiO.— 34 percent 


B.0;—28 

PbO—29 sy 
Al,O;— 7 
Na,O— 2 = 


1075 glass has satisfactory chemical and elec- 
trical properties. 

The alloy, referred to as 42 alloy, which matches 
this glass, has the composition : 


Iron— 57.8 percent 
Nickel— 41.5 
Manganese— 0.5 
Silicon— 0.2 
Carbon— <0.06. 





The thermal expansions of 42 alloy and 1075 
glass are shown in Fig. 10. It will be noted that 
the expansion in the range 0-300°C is essentially 
the same as that of fernico and 705-AO glass; 
but that the sealing point occurs at a lower 
temperature and much higher up on the curve, 
i.e., the expansion at the sealing temperature 
departs more from the low temperature expan- 
sion than in the case of 705-AOQO; the average 
expansion from room temperature to the sealing 
temperature being 1.4 times the 0-300°C value, 
as compared to 1.15 for the fernico—705-AO 
combination. This difference is characteristic of 
the two types of glass. As a result of this, the 
curves of glass and metal depart slightly in the 
range just below the sealing temperature, to join 
again at about 350°C. The stresses caused by 
this departure are not sufficient to cause breakage 
during cooling, and the seal is free from stress 
at all temperatures below 350°C. The photo- 
elastic stress pattern of a seal between 1075 
glass and 42 alloy, at room temperature, is 
shown in the insert in Fig. 10. 

A special problem encountered in making 
seals with 1075 glass is the appearance of lead 
at the interface between glass and metal, lead 
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oxide being reducible by ferrous metals. This can 
be avoided by coating the iron-nickel alloy with 
copper or platinum, which do not reduce lead 
oxide; or with a lead-free glass, for example, a 
pure borosilicate. We have obtained excellent 
results by “precoating” the metal with a thin 
layer of 705-AO glass. 


4. CHROME-IRON SEALS 


Seals between 26 percent chromium-iron alloy 
and Corning B-12 glass are highly strained. We 
have obtained satisfactory seals with this alloy 
and a glass of the following composition : 


Composition of Laboratory Glass No. 286 
SiO.—54 percent 
PbO—29 
K,0— 8 
Na,O— 5 
BaO— 4 





Figure 11 shows the relative thermal expansion 
of 26 chrome-iron and 286 glass. The insert 
shows the stress pattern of a seal between them. 


5. Iron SEALS 


Although the art of enameling is old, the 
glasses used as enamels are generally of lower 
expansion than the metal. This is desirable for an 
enamel, since it gives tangential compression of 
the glass, which prevents cracking, while the 
radial stress, which is tension, is small so long 
as the thickness of the enamel is small compared 
to the radius of curvature. But it is not desirable 
for a thick glass seal. The standard commercial 
glasses, also, are all much lower in expansion 
than iron, as can be seen from Fig. 9. 

For thick seals, as contrasted with thin 
enamels, a glass is needed whose average expan- 
sion between room temperature and the sealing 
point matches iron as closely as possible. A per- 
fect match, such as that between fernico and 
705-AO glass, cannot be had, since the expansion 
of iron is nearly linear. In addition to having 
the correct expansion, the glass must satisfy the 
usual requirements, for electrical applications, of 
high electric resistance and of chemical stability, 
especially with respect to water and acids. These 
requirements, as well as the expansion match, 
have been found to be satisfied by a glass of the 
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following composition : 


Composition of 542 Glass 
SiO.—45 percent 


K,O—14 
Na,O— 6 
PbO—30 
CaF2— 5 


Figure 12 shows the thermal expansion of this 
glass compared to iron; and the insert shows the 
stress pattern of a seal, annealed to 1°/min. 
It is seen that the seal is practically stress-free. 
Direct application of 542 glass to iron results 
in a layer of lead at the interface, as already 
noted in the case of 1075 glass and 42 alloy. 
This can be avoided by copper plating, or by 
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Fic. 12. Thermal expansions of 542 glass and iron, and 
photoelastic stress pattern of a seal between them. 


precoating with a leadless glass. We have had 
good results from precoating with standard soda- 
lime lamp glass. 


6. THE STRENGTH OF SEALS 


It is obvious that a seal which is entirely free 
from stress may be worthless, unless the glass 
adheres to the metal. This second half of the 
problem of sealing, viz., the question of adhesion, 
must be dismissed with only a brief discussion ; 
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not because it is less important, but because it 
has not yet been studied. A few rough measure- 
ments and some general observations are all that 
can be given at present. 

For the purpose of discussion, several tech- 
niques will be described: 


(1) Thin oxide bond 


The thickness of the oxide formed by heating 
the metal is carefully controlled: if too thin, it is 
dissolved by the glass; if too thick, it flakes off. 

An example is a combination of tungsten with 
Corning Nonex Glass. The tungsten wire is 
carefully oxidized in the flame. A thin tube of 
“wetting glass’ is then slipped over it and 
quickly sealed, thus preventing further oxidation. 
The coated wire may now be incorporated in any 
desired seal, care being taken not to injure the 
oxide layer by overheating. A certain “straw” 
color is a reliable index of correct oxide thickness. 

This may be termed the Ideal Sealing Tech- 
nique, since tungsten oxide dissolves in glass 
sufficiently to form a firm bond, but not so 
rapidly as to be completely absorbed; and the 
thin oxide layer may be expected to adhere 
firmly to the metal. The strength of the bond is 
not easily measured directly, since this sealing 
technique has been applied only to wires. But 
two indirect observations place upper and lower 
limits to the strength: (a) The bond is weaker 
than the glass, since tungsten-Nonex seals which 
are crushed in a vice break at the oxide-metal 
interface. (b) Tungsten-Pyrex seals, in which the 
radial tensile stress at the glass-metal surface is 
4.8 kg/mm?,® frequently do not break during the 
first few hours, but invariably break in less than 
a week. Hence the bond initially is stronger than 
4.8 kg/mm’. 


(2) Thick oxide bond 


The thick oxide layer behaves like a mis- 
matched enamel coat, which breaks loose from 
the metal during heating and cooling, because of 
the difference between its thermal expansion and 
that of the metal. Typical examples are seals to 
iron and copper, which over-oxidize easily. The 
danger may be minimized by careful sealing 


5 Reference 1, p. 407, Table VII. 
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Fic. 13. Form of seal used for testing the strength of 
bond between metal and glass. The studs, with the glass 
between them, are screwed into the jaws of a steel-testing 
machine. 
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technique, such as is employed in machine-made 
Dumet seals. Here two necessary conditions are 
fulfilled: (a) The wire is kept hot until covered 
with glass, so that cracking of the oxide layer, 
which is irreparable, does not occur before 
sealing. (b) The glass contracts more than the 
metal during cooling, and thus holds the oxide 
layer tightly to the metal, both during and after 
cooling. The strength of the bond is unimportant 
in this case, since the radial stress is compression. 


(3) Anchored oxide 


An example is 26 percent chromium-iron alloy, 
deeply etched by chemical pickling or by heating 
in wet hydrogen. The pores of the honeycomb- 
like surface support the oxide, which bonds 
firmly to the glass. Both lead glasses and lime 
glasses seal strongly to such a surface, and break 
in the glass rather than at the interface when 
crushed in a vice. Two direct measurements of 
the strength of bond between flat surfaces of 26 
chrome-iron and 286 glass gave tensile strengths 
of 4.4 and 6.8 kg/mm?. In contrast to this, 
unetched 26 chrome-iron gives a very weak 


bond. 


(4) Oxide-free seals 


Large seals usually require such a long heating 
time that the oxide layer is completely dissolved 
in the glass when the seal is finished. This is 
usually true of fernico-705-AO seals. Measure- 
ment of the strength of such seals will therefore 
be a test of the bond between glass and clean 
metal. 

We have made rough tests of the strength of 
bond in this and a few other sealing combinations 
by joining the heads of two tension bolts with a 
layer of glass, as shown in Fig. 13. The glass was 
kept smaller in diameter than the studs, or was 
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ground off flush with the metal if it protruded, 
so that the tension was purely normal. All seals 
were annealed and were essentially free from 
strain when tested. 

The strength was measured in a standard 
steel-testing machine. Tests of 28 seals between 
fernico and 705-AO glass, of 3-inch diameter, 
gave values ranging from 1.14 to 3.36 kg/mm*, 
the average being 1.48 kg/mm. All but 7 of 
these broke in the glass at the center, midway 
between the seals; in these cases the values 
obtained are those of the strength of the glass; 


they are not the bond strength, but only lower 
limits. For the 7 which broke, wholly or partially, 
at the glass-metal interface, the average strength 
was 2.23 kg/mm’. 

A few seals of this same type were made 
between 42 alloy, which had been precoated with 
705-AO glass, and 1075 glass. Tensile tests of 8 
of these seals gave strengths varying from 1.56 
to 3.05, the average being 2.27 kg/mm?. This is 
somewhat higher than the value found for the 
fernico-705-AO seals, and may indicate a cer- 
tain advantage in precoating. 





Diffusion Calculations 


Interrelation Between Two Solutions of the Fourier Equation 


HERMAN A. LIEBHAFSKY 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received May 1, 1941) 


An examination of the interrelation between the sine series solution (for a slab with sealed 
periphery) and the error function solution (for a semi-infinite body) of the fundamental Fourier 
equation has yielded the following results: (1) Distribution functions according to the two solu- 
tions begin to diverge soon after diffusion has begun. (2) The error function yields an excellent 
simple approximation, valid until half the material initially present has been lost, for the loss of 
material from a slab with sealed periphery. (3) A simple, but extremely good, approximation is 
given for an infinite exponential series in the region of slow convergence. (4) Two simplifications 
in the calculation of diffusion constants have been given. (5) While the discussion has been 
restricted to diffusion it is applicable, with obvious modifications, to other cases of flow governed 


by the fundamental Fourier equation. 


N the evaluation of measured losses of 

plasticizers from polyvinyl chloride plastics! 
it proved expedient to investigate the interrela- 
tion of two solutions for the fundamental Fourier 
equation 


dc/dt= Dd*c/dx* (1) 


(D is the diffusion constant, whose dimensions are 
Pt; units need not be specified) for the linear 
diffusion in the x direction of a substance at 
volume concentration c. As a result, several 
matters of more general application were dis- 
covered, and these will be presented here. 


1 An account of this investigation by H. A. Liebhafsky, 
A. L. Marshall and F..H. Verhoek will appear in the Journal 
of the American Chemical Society. 
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The two solutions are well known.? The first, 
C, ‘Co=2 nf 
0 


the ‘‘error function”’ solution, applies to a semi- 
infinite solid initially of uniform concentration Co, 
whose one plane face (situated at x=0) is 
brought at t=0 to c=0, and maintained in this 


z/2(Dt)} 


exp (—*)dg 
=erf (x/2(Dt)'), (2) 


2 Cf. W. E. Byerly, An Elementary Treatise on Fourier's 
Series, etc. (Ginn and Company, 1893), where Eq. (2) is given 
on p. 85 as the solution of example 1; and Eq. (3) as the 
solution of example 1, p. 106. For values of the “error 
function,” or probability integral (Eq. (2)), see B. O. Peirce, 
A Short Table of Integrals (Ginn and Company, 1910), 
second edition, p. 116. 


707 








condition thereafter. The second solution, 


c/co=4/xLe-” sin rx/a+ }fe-™” sin 34x/a 


+ te-*5¥ sin 5ax/a+---] (3) 
where 
y= Dr'*t/a*=rr*/4 (3a) 
applies to a slab with sealed periphery, a units 
thick, initially of uniform concentration Co, 
whose two faces are brought at t=0 to c=0, and 
maintained in this condition thereafter. 


COMPARISON OF DISTRIBUTION FUNCTIONS 


The distribution functions (variations of c¢ 
with x) are easily obtained from Eqs. (2) and (3), 
provided D, t, and a are known. Assuming D=1 
and a=2 makes ¢ equal to r (a variable intro- 
duced for convenience), and simplifies the com- 
parison. Before t=0, c=co whatever the value of 
x, so that the initial distributions are identical. 
They will remain sensibly identical thereafter so 
long as ¢ remains sensibly equal to cy at x= 1, the 
midpoint of the slab. At t=0.2=7, they diverge 
appreciably, as Fig. 1 shows. (Here the portion 
from x=0 to x=1 of the semi-infinite body 
represents also an outer half of the slab; the 
distribution function for the slab is, of course, 
symmetrical about the center.) The semi-infinite 
body has the higher c/co (0.886 against 0.772 at 
x =1) because of material that has flowed in from 
the semi-infinite reservoir beyond x=1. The 
error function solution (Eq. (2)) cannot be used 
to calculate accurately the distribution within a 
slab very long after diffusion has begun. 


COMPARISON OF AMOUNTS OF MATERIAL LOST 


The amount of material AQ lost from unit area 
in the interval t=0 to t=¢ can be calculated from 
Eqs. (2) or (3) after integrating with respect to ¢ 
or x. For the semi-infinite body, the former 
integration is more convenient: 


t 
Q.-0.=a0= f D(dc/dx) -- pdt =co2n'(Dt)!. (4) 


For the slab, integration with respect to x 
defines an average concentration @ (the ratio R is 
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Fic. 1. Distribution function at r =0.2 according to Eq. (2) 
(upper curve) and Eq. (3) (lower curve). 


introduced for convenience): 


é/co= R=1 af (c/co)dx 


= 8/m*[e-¥+ 146-94 + 14 5e774 +--+]. (5) 
Since, for a slab of volume V, 
Qo=coV; Qi=eV; 
R=Q:/Qo=1—(A4Q/Qo), 


R values can be calculated also from Eq. (4). 
This calculation is illogical to the extent that é 
for the entire semi-infinite body will never be 
appreciably less than co; what the calculation 
gives (approximately) is R for a unit cube that is 
part of the semi-infinite body, with one face at 
x=0 (Fig. 1). We shall call this calculation of R 
the “error function approximation” for the slab. 
For this unit volume, 


R=1-—AQ/co 
= 1—(2/2'(Dt)')(approximately). (7) 


6) 


Values of R from Eq. (5), tediously computed by 
taking enough terms for five significant figures, 
and from the much simpler Eq. (7), are given in 
Table I. 
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Calculational errors are obviously responsible 
for the unsystematic difference of 1 in the fifth 
significant figures of columns 2 and 3, Table I, 
for r=0.12 and less. (Column 4 will be explained 
later.) This remarkable agreement is the more 
surprising in view of the divergence of the 
distribution functions for the two solutions. Even 
for r=0.2, where Fig. 1 shows this divergence to 
be marked near x=1, the R values differ only by 
1 part in a thousand, which suffices for most 
diffusion measurements. (In the work on plastics, 
weight changes to 1 part in 500 of the plasticizer 
could be detected.) In other words, the simple 
Eq. (7) can replace the cumbersome Eq. (5) for 
all ordinary diffusion measurements on a slab 
with sealed periphery so long as not more than 
half the material originally present is lost. 

Equation (7) is thusan excellent approximation 
for the exponential series, Eq. (5), at values of y 
below 0.3, where evaluation of the series is 
relatively tedious. Furthermore, this mathe- 
matically interesting approximation is somewhat 
unlikely of discovery by mathematical methods 
alone. 
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TABLE I. Calculated values of R for the slab. 








ERROR FUNCTION 








Exact SOLUTION APPROXIMATION ERROR FUNCTION 
T (Eo. (5)) (Ec. (7)) SOLUTION 
0.01 0.88717 0.88716 0.8872 
02 .84041 .84042 .8404 
04 .77432 .77432 .7744 
.06 .72360 .72361 .7240 
.08 .68084 .68085 .6824 
.10 .64318 .64318 .6470 
R .60913 .60912 -- 
16 54876 54865 5814 
20 49591 49537 5554 


The reason for this astonishing agreement is 
clear from Fig. 1. The area I represents material 
that has flowed from the semi-infinite reservoir 
into the unit cube under consideration; the area 
II represents the excess of material resulting from 
the divergence of the two distribution functions 
(Eqs. (2) and (3)). So long as these two areas are 
sensibly equal, the concentration gradients for 
the two solutions will also be sensibly equal near 
x=0. The amount of material that has left the 
cube through the face at x =0 will consequently 
be identical for the semi-infinite body and for the 
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slab; hence the agreement of columns 2 and 3, 
Table I. 

Under the conditions of Fig. 1, the true R for 
the slab (Eq. (5)) is 1 — (area IV+-area II); R from 
Eq. (7) is 1—(area I1V+area I); the R actually 
realized for the unit cube attached to the semi- 
infinite body is 1—(area IV), and it can be ob- 
tained by adding values of area I to the proper 
values of R from Eq. (7). Area I has been graphi- 
cally evaluated for a number of rt values; column 
4, Table I, and the upper curve, Fig. 2, give the 
corresponding exact values of R for a unit cube 
attached to the semi-infinite body. These R 
values are of interest mainly because they show 
that an exact calculation of R from the error 
function is useless as an approximation for the 
case of the slab. 

THE CALCULATION OF DIFFUSION CONSTANTS 

The use of the error function approximation 
facilitates the calculation of diffusion constants 
for the slab series (Eq. (5)) and, by way of 
reciprocation, this series can be used to simplify 
the calculation of diffusion constants from the 
error function solution. 

To evaluate D from R values in agreement 
with Eq. (5), Sherwood* has used a special 
ordinate scale that gives a straight line when such 
values are plotted against ¢ or 7 as abscissa. Our 
experimental results required the computation 
of a scale considerably more detailed than 
Sherwood’s.* Equation (5) converges very slowly 


3 (a) T. K. Sherwood, Ind. Eng. Chem. 21, 12 (1929). (b) 
Cf. A. B. Newman, Trans. Am. Inst. Chem. Eng. 27, 310 
(1931). 

‘Figure 2 was prepared as follows. r=0.6, the highest 
abscissa of practical interest, was laid off 30 cm from the 
origin on millimeter graph paper. R=1.000, the ordinate 
for t=0=7, was laid off 24 cm from the origin, and a 
straight line, the reference line, was drawn through the 
two points thus plotted. The ordinate scale that will make 
this reference line satisfy Eq. (5) for D=1, a=2, is now 
required. It was obtained by solving Eq. (5) for values of 
rt from 0.60 to 0.02 in steps of 0.02, and for the value 
7=0.01. The resulting R values are given in Fig. 2, in 
which a change of 0.02 in 7 is a change of 10 mm on the 


at small 7; at r=0.01, six terms must be calcu- 
lated to establish R to three significant figures. 
The use of Eq. (7) resolves this difficulty for R 
values above 0.5; between 0.5 and 0.45, only two 
terms of Eq. (5) are required for three significant 
figures in R; thereafter, only one term is neces- 
sary. For most measurements, therefore, it 
suffices to use Eq. (7) until half the material 
initially present in the slab has been lost, and the 
first term of Eq. (5) thenceforward. A comparison 
of the three types of R values in Table I is given 
at the higher r’s by Fig. 2; the error function 
solution curve tends to be horizontal when the 
rates of flow across the boundaries x =0 and x=1 
are roughly equal. 

A simplification follows of the calculation of D 
from measurements that fit the error function 
solution, or from R values near unity for the slab. 
The usual plot of weight losses (AQ’s) against ft} is 
made. The slope of the resulting straight line 
(Eq. (4)) divided by Qo defines k = (1— R) /t*. But 
Fig. 2 relates R to r; at r=0.5, for example, 
1— R=0.764. Now, by Eq. (3a), 


(0.764/k)?=t=a?/4X4X1/D, 
D=(ak)?/4.670. (8) 


This combination of the sine series and error 
function solutions of the Fourier equation has at 
least two advantages: R is dimensionless so that 
concentration units disappear from the calcula- 
tion of D; and the areas of the specimens need 
not be known, provided no appreciable dimension 
changes occur during the experiments. 

I am indebted to my colleagues in the General 
Electric Company, especially to Drs. Hillel 
Poritsky, R. P. Johnson, and Lewi Tonks, for 
valuable discussions bearing on the foregoing 
material. 


original abscissa scale, and corresponds to a change of 
8 mm on the original ordinate. 

If experimental R values when plotted against ¢ lie on 
a straight line different in slope from the reference line, 
Eq. (5) is obeyed; but the conditions D=1, a=2, are not 
both satisfied. D at any such R can be calculated by Eq. 
(3a), ¢ being the experimental abscissa for this value of R; 
7, that from the reference line. 
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